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1.  INTRODUCTION 


In  recent  years,  considerable  effort  has  been  expended  at  The  Aerospace 
Corporation  Aerophysics  Laboratory  to  design  theoretical  models  for  computer 
calculations  of  the  behavior  and  performance  of  HF  chemical  lasers,  both 
pulsed  and  continuous.  In  order  to  perform  such  calculations,  the  rate  coef- 
ficients of  the  varioucr  chemical  reactions  and  energy  transfer  processes  that 
occur  must  be  known,  or  at  least  estimated.  To  this  end,  the  Chemical 
Kinetics  Department  has  attempted  to  maintain  an  up-to-date  compilation  of 
the  rate  coefficients  thet  reflects  the  best  current  information  available  from 
published  and  unpublished  theoretical  and  experimental  work.  When  direct 
information  has  not  been  available,  theoretical  techniques  for  estimating  rate 
coefficients  ba3ed  on  various  a priori  principles  or  estimates  based  on  com- 
parison w.th  related  reactions  with  known  rate  coefficients  arc  applied. 

This  report  is  the  sixtnof  a series  prepared  by  the  Chemical  Kinetics 
Department  during  the  past  four  years  in  which  the  kinetics  of  hydrogen  halide 
laser  systems  are  reviewed  [Cohen  (1971,  1972,  i974);  Cohen  and  Bott  (1975, 
1976)].  For  the  users'  convenience,  each  report  is  self-contained,  obviating 
recourse  to  the  preceding  reports.  However,  as  additional  information  has 
become  available,  or  as  interests  have  shifted  slightly,  certain  issues  dealt 
with  in  earlier  reports  did  not  need  to  be  covered  again,  and,  to  this  extent, 
these  reports  are  not  completely  cumulative  in  scope. 

Hydrogen-halide  kinetics  and  energy  transfer  is  a field  in  which  there 
is  great  activity,  and  new  results  are  continually  becoming  available.  There- 
fore, it  is  inevitable  that  information  in  this  report  will  soon  be  out-of-date, 
perhaps  before  it  is  disseminated.  To  minimize  this  problem,  reference  is 
made  to  unpublished  data  and  reports  when  appropriate.  The  assistance  of 
many  colleagues  who  have  provided  preprints  cf  their  work  is  gratefully 
acknowledged.  The  literature  survey  for  this  review  was  completed  in 
June  1975. 
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The  scope  of  this  survey  includes  all  reactions  of  interest  in  an 
H^-F^  chemical  laser  system.  Thus,  chemical  reactions  and  energy-trans- 
fer processes  involving  other  hydrogen-  or  fluorine-containing  compounds 
have  been  omitted.  In  most  current  chemical  kinetics  computer  programs, 
rate  coefficients  are  entered  in  the  form  ATncxp( -E/RT).  In  this  ~eport, 
therefore,  all  coefficients  have  been  made  to  fit  this  form,  even  when  data 
or  theory  suggest  a different  temperature  dependence.  The  primary  purpose 
here  is  to  present  data  for  computational  purposes,  rather  than  to  unravel 
the  theoretical  significance  of  the  findings.  Joint  Army,  Navy,  and  Air 
Force  (JANAF)  data  have  been  used  for  thermochemical  quantities,  with 
Ka  k = k^/k^  uSed  to  express  the  equilibrium  conntant  in  terms  of  the  rate 
coefficients.  Units  of  cubic  centimeters,  moles,  seconds,  and  calories  are 
generally  used  throughout.  In  fitting  data  with  analytic  expressions,  particu- 
lar attention  has  been  given  to  the  temperature  range  300  to  1000  K;  higher 
temperatures  are  not  of  practical  interest,  and  lower  temperatures  are 
largely  devoid  of  data.  Extrapolation  of  the  recommended  rate  coefficients 
outside  of  this  range  should  be  carried  out  with  caution. 
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H.  DISSOCIATION-RECOMBINATION  REACTIONS 


A.  11  DISSOCIATION 

Hydrogen  dissociation  and  recombination  rates  have  been  widely 
studied  since  the  first  room-temperature  measurements  were  made  in 
1929.  The  data  have  been  summarized  and  reviewed  critically  by 
Baulch  et  al.  (1972).  Because  these  reactions  are  only  of  slight  importance 
in  laser  systems,  which  are  generally  at  low  (<500  K)  temperatures,  all 
the  experimental  evidence  is  not  tabulated  here.  Shock-tube  studies  of  the 
following  reactions  have  been  made  by  Jacobs,  Giedt,  and  Cohen  (1967) 
in  the  temperature  range  2900  to  4700  K. 

H + H + Ar  it  H2  + Ar  (la) 

H + H + H-  = Hn  + H0  (lb) 

Luc, 

H + H + H^H2  + H (lc) 

Their  review  of  previous  shock-tube  data  showed  that  the  results  of  various 

experimenters  disagreed  by  factors  of  about  3,  5,  and  10,  respectively,  for 

Ar,  H,  and  as  chaperone  gases.  (The  shock-tube  experiments  actually 

involved  the  measurement  of  hydrogen  dissociation.  However,  it  has  been 

customary  to  report  the  results  in  terms  of  a recombination  rate.)  V all 

three  cases,  their  results  lie  in  the  middle  of  the  range  of  values.  They 
18  - 1 

reported  k.  -10  T , k 4 , 2.5k  , and  k,  = 20  k . . Their  values 

r la  lb  la  lc  la 

for  and  k^,  extrapolated  to  room  temperature,  yield  3.3  and  8.  2 X 101 
cc/mol-sec,  respectively.  In  four  decades,  the  room  temperature  value 
of  k^  has  ranged  from  0.2  to  72  X 10*^.  Most  recently,  Ham,  Trainor, 
and  Kaufman  (1970)  [see  also  Trainor,  Ham,  and  Kaufman  (1973)]  have 
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measured  k^  and  k^  in  a flow  system  using  an  isothermal  catalytic  probe 
over  the  temperature  range  of  77  to  300  K and  found  temperature  dependen- 
cies of  t~0*81  and  approximately  T’0*^1,  respectively,  with  values  at 

1 S 2 2 

300  K of  3.3  and  2,9  X 10  cc  /mol  -sec.  Their  results,  extrapolated  to 

the  temperature  range  of  the  experiments  of  Jacobs,  Giedt,  and  Cohen  (1967)f 

agree  with  the  results  of  the  latter  workers  within  35%,  Because  tne  lower 

temperatures  are  of  more  interest  in  laser  chemistry,  the  results  of  Ham, 

Trainor,  and  Kaufman  (1970)  are  weighted  heavily.  The  expressions 

k.  = 6.2  X 10 17  T"0, 95  and  k..  = 9.4  X 10l6T’°’61  fit  both  high-  and  low- 
la  lb 

temperature  data  within  about  10%.  The  numerous  results  are  not  tabulated 
here;  further  citations  are  given  by  Baulch  (1972). 

Of  the  several  room- temperature  measurements  of  k^,  the  most 

reliable  seems  to  be  that  of  Bennett  and  Blackmore  ( 1968).  They  obtained 

15  2 2 

an  upper  limit  of  2,5  X 10  cc^/mol  -sec  in  a discharge  flow  system  using 
esr  for  atom  detection,  suggesting  that  kjc  has  a maximum  value  some- 
where in  the  temperature  range  1000  to  3000  K and  is  relatively  unimportant 

14  1/2 

at  low  temperatures.  The  expression  k ^ = 1.2X10  T is  consistent 
with  the  upper  limit  of  Bennett  and  Blackmore  (1968)  and  agrees  with  the 
low-temperature  end  of  the  shock-tube  experiments  of  Jacobs,  Giedt,  and 
Cohen  (1967).  Therefore,  it  is  probably  a reasonable  expression  in  the 
temperature  range  of  interest*  although  it  is  not  useful  above  3000  K.  This 
expression  is  recommended  in  Appendix  1,  subject  to  modification  when  more 
experimental  data  become  available. 

B.  F2  DISSOCIATION 

A thorough  critical  review  of  the  dissociation-recombination  of 
and  Cl0  was  published  by  Lloyd  (1971);  repetition  of  his  work  has  been  kept 
to  a minimum  here. 

In  recent  years,  several  shock-tube  determinations  have  been  made 
of  the  rate  of  dissociation  (Table  1). 

F>  + M«F+F  + M (2) 
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Althougl  the  numerical  values  of  the  rate  coefficient  have  generally  been 
in  agreement,  within  the  experimental  uncertainties  of  the  various  deter- 
minations, the  temperature  dependencies  of  the  rate  coefficient  have  been 
noticeably  inconsistent.  The  various  results  are  summarized  in  Table  1. 
The  rate  coefficients  calculated  from  shock-tube  data  are  plotted  (Fig.  1); 
the  recombination  rate  coefficients  have  been  calculated  and  extrapolated  to 
lower  temperatures  (Fig.  2).  The  extent  to  which  the  F^  bond  energy 
exceeds  the  average  experimental  activation  energy  is  disconcerting.  This 
discrepancy  becomes  important  when  calculating  the  rate  coefficient  for 
F + F recombination  from  ^ and  and  extrapolating  the  results  to  low 

temperatures.  The  problem  is  compounded  by  the  uncertainty  in  the  bend 
dissociation  energy  of  F->. 


Although  a critical  review  of  the  bond  strength  of  F^  is  outside  the 
scope  of  this  report,  it  is  evident  that  the  possible  uncertainty  in  such  a 
fundamental  datum  leads  to  ambiguities  in  several  thermochcmical  and 
kinetic  values  important  to  this  study.  In  1968,  a National  Bureau  of 
Standards  committee  recommended  37.  76  keal/mol  foi  D298^2^  since 
1965,  the  JANAF  committee  has  recommended  37.72  ±0.8  kcal/mol.  (The 
dissociation  energy  at  0 K,  Dq,  in  either  case  is  1.0  kcal/mol  smaller.  ) 
Here,  38  kcal/rrol  is  used;  in  view  of  the  uncertainties,  a third  significant 
figure  is  unnecessary. 


With  a value  of  38  kcal/mol  for  D^gfF^),  all  of  the  shock-tube  values 

for  the  activation  energy  of  F^  dissociation  in  an  Ar  bath  lead  to  a strongly 

negative  temperature  dependence  for  F-atom  recombination,  which  results 

in  unreasonably  large  values  near  300  K.  Therefore,  for  calculations,  we 

prefer  the  theoretical  value  calculated  by  Shui,  Appleton,  and  Keck  (197  1). 

Their  results  agree  numerically  with  the  shock-tube  data  of  Johnson  and 

Britton  (1964)  but  lead  to  rn  activation  energy  of  35.  1 kcal/mol  and  a pre- 

| > 

exponential  factor  of  10  . Until  direct  measurements  near  300  K are 

^ r 13  7 

available,  10  * exp(-35,  100/RT)  is  recommended  which  implies 
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Ar  1 4 0 

that  k . at  300  K is  iO  . A recent  report  by  Ganguli  and  Kaufman  (1974) 
suggests  log  k ^ **  13.  5 at  300  K,  which  is  smaller  than  the  recommended 
value  for  k ^ by  a factor  of  S3. 

As  a result  of  their  recent  work,  Breshears  and  Bird  (1973)  have 
concluded  that  k^/k^r  is  approximately  2.7  in  the  temperature  range  1400 
to  2600  K.  However,  an  examination  of  Fig.  1,  in  which  alt  the  shock-tube 
dissociation  data  are  plotted,  shows  that  their  data  are  not  greatly  sensitive  to 
this  ratio;  therefore,  considerable  uncertainty  must  be  assigned  to  this  value. 
Data  on  the  rate  of  I ^ dissociation  with  M = H,  H?,  or  F are  not  available. 
From  an  analogy  with  the  findings  incases  of  other  homonuclcar  diatomics  such 
as  Cl-,,  Br?.  H?,  0?,  and  N-,,  it  is  expected  that  F atoms  will  exhibit  a large 

L l l p<: 

efficiency  relative  to  Ar;  k^  = 10  k^  is  assumed  in  this  report. 

If  it  is  assumed  that  the  primary  recombination  mechanism  for  F 
atoms  is,  as  is  believed  to  be  the  case  for  I atoms,  the  radical-molecule 
mechanism 


F + M 

FM*  + F — F2  + M 

then  some  qualitative  statements  can  be  made  regarding  the  relative  effi- 
ciencies of  H and  HF  as  chaperones.  For  example,  H atoms  would  be 
expected  to  be  somewhat  more  efficient  than  Ar  because  the  strong  H-F 
bond  would  favor  formation  of  FM*,  where  M = H.  The  effect  on  HF 
chaperon*  efficiency  would  be  much  less  pronounced.  Therefore,  in  this 
report,  it  is  tentatively  assumed  that  H atoms  are  three  times  as  efficient 
as  Ar  and  that  HF  is  equal  to  Ar  in  efficiency.  Other  species  are  assumed 
to  have  the  same  efficiency  as  Ar. 

Figure  2 is  significant  in  demonstrating  that  a relatively  small  dis- 
crepancy in  k at  shock-tube  temperatures  can  lead  to  an  unacceptably  large 
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divergence  of  extrapolated  values  at  room  temperature.  When  the 
uncertainty  in  the  F^  bond  dissociation  energy  (and,  therefore,  the  equilib^ 
rium  constant)  is  taken  into  account,  the  necessity  for  a reliable  direct 
measurement  at  300  K becomes  apparent.  For  the  present,  however,  it 
must  be  accepted  that  the  room- temperature  recombination  rate  is  uncer- 
tain by  at  least  an  order  of  magnitude. 

C.  Hr  DISSOCIATION 

The  rate  of  HF  dissociation 

HF  + M — H + F + M (3) 

was  first  measured  by  Jacobs,  Giedt,  and  Cohen  ( 1 96 S > in  the  shock  tube  at 

temperatures  of  3800  to  5300  K.  They  reported  a rate  coefficient  for  Ar  as 

, . . f . Ar  1A19. 05^-1  , ^ /rv~.  , Ar  IA22.7l 

chaperone  in  two  farms:  k-  = 10  T cxp(-D  /RT)  or  k-  = 10 

2 ~ o j 

T exp(-DQ/RT).  where  Dq.  the  dissociation  energy  at  0 K from  v = 0,  is 

currently  taken  as  135.  1 kcal/mol.  (At  the  time  of  their  study,  a slightly 

smaller  value  of  134.  1 kcal/mol  was  accepted.  ) Subsequently,  Blauer  (1968) 

performed  similar  shock-tube  experiments  over  the  temperature  range 

3700  to  6100  K and  obtained  a rate  coefficient  smaller  by  about  a factor  of  2, 

with  a pre -exponential  factor  of  10’  ^’^T  The  two  studies  were  reviewed 

by  Brown  (NBS  Report  No.  10-635),  who  concluded  that  there  was  no  basis 

for  preferring  one  result  to  the  other.  Brown  assembled  all  the  data  of 

both  studies  in  one  plot  (Fig.  3). 

Because  of  advances  in  the  knowledge  of  related  rate  coefficients,  it 
is  worthwhile  to  re-examine  the  data  of  Jacobs,  Giedt,  and  Cohen  (1967)  to 
determine  if  their  conclusions  are  still  valid.  They  deduced  the  value  of  by 
computer  modeling  the  profile  of  (HF)/(HF)  f where  the  denominator  is  the 
initial  HF  concentration.  Although  their  model  included  all  five  reactions 
occurring  in  the  system,  i.e.,  HF  dissociation- recombination , 
dissociation- recombination , F^  dissoc iation- recombination , H + HF,  and 
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Fig.  3.  Experimental  Measurements  of  HF  Dissociation  Rate 
Coefficient  Versus  IO'Vt 
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F + HF,  they  found  that  two  of  the  five  reactions  did  not  contribute 
measurably  to  the  overall  kinetics  and  could  be  ignored.  The  system  of 
significant  reactions,  therefore,  consists  of 


HF  + M.  ,^H+F  + M.  . 
3,i  3,i 


(3-i) 


H + HF=H2+  f 


(4) 


H~  + M4  .5SH  + H + M.  . 

2 1,1  l,i 


(1-i) 


Reaction  (1)  was  considered  as  three  separate  processes,  (la),  (lb),  and 

(lc),as  defined  earlier,  where  a,  b,  anJ  c represent,  respectively,  M = Ar, 

H^,  and  H.  was  assumed  to  be  a weighted  average  over  all  species 

present  because  no  experimental  basis  was  found  for  assigning  different 

efficiencies  to  the  different  species  (principally  Ar  and  HF).  The  value 

1 3 

they  used  for  was  10  exp( -35000/RT).  which  with  currently  accepted 

thermochemical  values,  implies  a value  for  the  reverse  rate  coefficient, 

14  «n  471 

k 4,  of  5.5  X 10  T * exp{-6490/RT).  At  presenc,the  most  probable 


-4' 


value  for  k ^ is  1.6  X 10^  exp(- 1600/R  T)  (Section  II),  wnich  yields  a value 

at  4000  K that  is  27.5  times  larger  than  that  used  by  Jacobs,  Giedt,  and 

5 36 

Cohen  (1967),  and  at  300  K,  it  is  10  * times  as  large. 

For  kj^,  Jacobs,  Giedt,  and  Cohen  (1967)  used  10*^*^T~*,  which  is 

about  twice  as  large  as  the  current  best  estimate,  as  discussed  earlier. 

Furthermore,  since  they  scaled  k..  and  k4  relative  to  k,  (they  used 

7 lb  1c  la  1 7 

k^  = 2.47  kja#  kjc  = lOkj^),  these  rate  coefficients  are  also  subject  to 

revision.  Current  best  estimates  put  at  approximately  2.5  k^a  and 

klc  — 20  kj^.  Thus,  their  kj^  should  be  decreased  by  a factor  of  2,  but 

their  k4  is  correct.  When  these  revisions  were  made  in  the  values  for  k. 

lc  1 

and  k^  and  the  computer  profiles  rerun,  the  agreement  with  the  experimental 
profiles  was  not  quite  as  good.  Several  variations  in  k^  were  tried,  but  none 
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was  found  that  gave  as  good  agreement  with  the  data  throughout  the 

temperature  range  as  did  the  earlier  set  of  rate  coefficients.  However,  the 

best  overall  fit  is  obtained  with  kl* ^ = 8 X 10^T  * exp(-  111, 000 /RT)  cc/mol 

* Hp 

sec#  which  is  about  30%  larger  than  the  previous  value  of  at  4000  K and 
about  30%  smaller  at  5000  K.  However,  because  of  the  larger  deviations 
between  this  revised  rate  coefficient  and  the  experimental  data,  there  is 
some  doubt  that  this  expression  is  more  useful  than  the  previous  one  when 
extrapolated  into  the  regime  of  lower  temperatures,  which  is  of  greater 
interest.  The  large  discrepancy  between  this  revised  apparent  activation 
energy  of  1 1 1 keal/mol  and  the  HF  bond  energy  will  lead  to  unreasonably 
large  dissociation  rates  near  room  temperature.  (This  problem  was  dis- 
cussed in  connection  with  F ^ dissociation  in  Section  I.  D.  ) Therefore,  for 
the  present,  the  earlier  value  obtained  by  Jacobs,  Giedt,  and  Cohen  (1967) 
is  recommended  because  of  its  greater  use  at  low  temperatures. 


111.  PUMPING  REACTIONS 


A.  F + H2—  HF  + H 

1.  OVERALL  REACTION  RATE 

Until  recently,  there  have  been  few  reports  of  direct  experimental 
measurements  of  F-atorrv  reactions,  partly  because  of  the  handling  difficulties, 
and  partly  because  of  the  complications  resulting  from  the  presence  of  molec- 
ular in  most  cases.  Experimental  studies  of  the  F ► reaction  fall  into 
three  categories:  (a)  direct,  absolute  rate  measurements;  (b)  relative,  com- 

petitive fluorination  measurements;  and  (c)  complex  systems  in  which  the  rate 
coefficient  is  inferred  by  detail  modeling. 

a.  Direct,  Absolute  Measurements 

Five  sepai  ate  studies  fall  into  this  category.  In  the  earliest,  Homann 

et  al.  (1970)  produced  F atoms  in  a flow  system  by  the  reaction  N ^NF^^N^  + 2F 

and  then  introduced  downstream.  Ir.  some  experiments,  the  disappearance 

of  F atoms  was  monitored  by  titration  with  C1NO;  in  others,  the  appearance  of 

HF  was  followed  with  mass  spectrometry.  From  measurements  at  400,  J60, 

14 

and  300  K,  an  Arrhenius  expression  of  k.  = 1.6  x 10  cxp(  - 1600/RT)  cc/mol- 

* 13 

sec  was  obtained  with  a room-temperature  value  of  1.  1 x 10  . There  arc 

a few  questionable  aspects  to  this  procedure.  The  possible  disadvantages  of 
the  CINO  titration  technique  have  been  discussed  by  Clyne.  McKcnney,  and 
Walker  (1973).  Another  question  concerns  the  effect  of  excess  N atoms  of 
NF2  radicals  if  the  titration  reaction  producing  F atoms  is  not  carried  out 
under  stoichiometric  conditions,  if  NF^  is  present  in  excess,  it  is  possible 
that  H atoms  produced  in  Reaction  (4)  could  react  with  NF0.  increasing  the 
yield  of  HF.  Since  this  process  is  probably  fast  compared  to  React  ion  (4), 
the  calculated  k^  would  then  be  too  large  by  a factor  of  up  to  2.  (The  pos- 
sibility that  N + NF2  reacts  to  give  2NF,  which,  in  turn,  produce  N2  + 2F. 
should  not  be  discounted,  because  it  would  lead  to  kinetic  complications  if 
the  reaction  between  N and  NF2  is  not  carried  out  stoichiomctr ically . ) 


iTtSCEDlJG  BLANK- ItCT  FILL'D 


Dodonov  et  al.  (1971)  studied  the  reaction  in  a fast  flow  of  He,  with  F 

atoms  produced  by  a high-frequency  discharge  through  F^ , diluted  with 

He,  was  added  downstream,  and  the  reaction  system  was  probed  by  a mass 

spectrometer.  They  observed  nearly  complete  disappearance  of  F^  when 

the  discharge  was  turned  on.  with  roughly  comparable  concentrations  of 

V atoms  and  SiF^  (from  reaction  with  the  glass  walls)  being  formed.  If  the 

reaction  rate  of  H atoms  with  SiF^  is  negligible  at  room  temperature,  which 

It  probably  is,  then  these  experiments  give  a good  measure  of  the  value  of  k« , 

13  ** 

for  which  Dodonov  et  al.  (1971)  obtained  1.8  X iO  at  293  K , with  an  est i mated 
uncertainty  of  33%,  largely  resulting  from  uncertainty  in  the  F atom  concentra- 
tion measurement. 

Clyne,  McKcnncy,  and  Walker  (1973)  produced  F atoms  in  a microwave 
discharge  through  flowing  CF^  and  11c  and  reacted  the  mixture  with  H^.  The 
F-atom  concentration  was  monitored  by  titration  with  C^;  the  concentra- 
tion was  followed  by  mass  spectrometry.  The  rate  coefficient  at  298  K was 
1 3 

1.5  X 10  cc/mol-scc,  with  50%  uncertainty  because  of  considerable  experi- 
mental scatter.  In  the  same  system,  they  measured  tb^  rate  of  F + CH^, 

which  is  of  interest  in  connection  with  the  competitive  F + measurements 

1 3 ^ 

to  be  discussed  later,  and  obtained  3,6  X 10  , with  even  larger  scatter, 

Kompa  and  Wanner  (1972)  produced  F atoms  by  the  flash  photolysis  of 
WF^  at  298  K in  the  presence  of  in  a stationary  laser  cavity.  The  rcactior 
rate  coefficient  was  deduced  by  monitoring  the  decay  of  chemiluminescence 
from  vibrationally  excited  HF.  This  technique  requires  the  minimization  of 
the  vibrational-relaxation  process  and  of  the  secondary  reaction  of  H or 
with  WFx,  the  decomposition  products  from  WF^.  The  first  of  these  require- 
ments can  be  met  by  operating  at  sufficiently  low  pressures;  the  limitation  on 

pressure  then  becomes  the  sensitivity  of  the  optical  detection  system.  The 

1 3 

rate  coefficient  obtained  was  3.8  x 10  cc/mol-scc.  In  the  same  system, 
the  rate  coefficient  for  F + CH^  was  measured  to  be  4.3  X 10*  . A second 
possible  complication  is  not  so  easily  disposed  of.  Kompa  and  Wanner  (1972) 
operated  at  reactant  pressures  of  0.25  Torr  each  of  WF^  and  (or  CH^). 
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With  a 0.  b^o  photolysis  assumed,  the  H atoms  produced  in  Reaction  (4)  would 

9 

react  with  WF^  rather  than  combine  by  a three-body  process  if  — 10 

cc/mol-sec,  i.e.  f il  the  activation  energy  is  less  than  approximately  7 kcal/ 
mol.  For  reaction  with  WF^  the  rate  coefficient  would  have  to  be  approx- 
imately 100  times  larger  or  more,  since  (WF^J/JWF^  ] < 0. 0 1 . If  either  of 
these  processes  is  significant,  then  the  calculated  rate  coefficient  for  the 
F + reaction  could  be  too  large  by  as  much  as  a factor  of  2,  which  is  very 
possible  although  the  reactions  of  H with  WF^  and  WF^  have  not  been  studied. 
Also.  Igoshin,  Kulakov,  and  Nikitin  (1974)  have  pointed  out  the  possibility 
that  some  F atoms  in  the  above  work  were  t ranslationally  hot  because  of  the 
low  total  pressures,  also  resulting  in  a spuriously  high  value  for  k^.  Infor- 
mation is  not  available  to  evaluate  quantitatively  the  extent  to  which  this  might 
occur.  However,  a comparison  of  the  average  bond  strength  in  WF^  ( 1 Z t kcal/ 
mol)  with  the  peak  in  the  WF^  absorption  spectrum  Incar  1800  A,  according 
to  Kompa  and  Wanner  (1972)]  suggests  that  the  F atoms  formed  in  the  photo- 
lytic  act  could  easily  have  30  keal/mol  excess  translational  energy,  which  they 
cannot  lose  ueforc  collision  with  (since  collision  with  WF^  will  not  remove 
significant  translational  energy).  This  excess  of  translational  energy  could  be 
the  most  serious  flaw  in  the  experiment. 

Hon,  Axworthy,  and  Schneider  (1973)  studied  the  reaction  of  F atoms 
with  a variety  of  potential  laser  system  fuels,  including  and  CH^.  Rate 
coefficients  were  obtained  in  a flow  system  with  F atoms  produced  by  thermal 
dissociation  of  F^  and  monitored  by  csr.  The  rate  of  reaction  of  F with 
at  298  K was  found  to  be  1 X 10**  cc/mol-s^c;  that  with  CH^  was  2,6  X lO1^. 

The  five  determinations,  thus,  range  from  1 to  4. 3 X 10**  cc/mol-scc 
and  arc  listed  in  Table  2.  At  present,  the  value  of  Clyne  ( 1973)  seems  to  be 
most  reliable,  in  spite  of  the  large  uncertainty  the  authors  place  on  their 
result.  Thus,  a value  of  = 1.5  X 10**  cc/mol-scc  is  suggested  by  the 
four  direct  absolute  measurements. 

b . Competitive  Studies 

The  first  relative  measurements  were  made  by  Mercer  and  Pritchard 
(1959).  In  a thermal  reactor  over  the  temperature  range  298  to  423  K,  they 
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Table  2.  Absolute  Measurements  of  F + H HF  + H 
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studied  the  competition  between  and  Cll^  for  F atoms  initiated  by  the 
photolysis  of  Fp.  The  consumption  of  reactants  was  monitored  by  mass 
spectrometry.  The  results,  together  with  the  other  competitive  measure- 
ments, are  shown  in  Table  3.  The  results  of  competitive  studies  in  Fig.  4 
were  calculated  with  the  value  for  F + Cll^  recommended  by  Fettis  and  Knox 
( 1964)  assumed. 


Kapralova,  Margolin,  and  Chaikin  (1970)  studied  the  competitive  fluori- 
nation  between  CH^  and  and  between  H^  and  Dp  in  spherical  heated  flasks, 
diluted  with  He;  the  products  HF  and  DF  were  analyzed  by  csr.  The  studies 
were  conducted  at  77,  173,  and  293  K for  CH^  and  Dp  and  223,  253,  293,  and 
353  K for  Hp  and  Dp.  The  Dp/CH^  data  arc  reasonable,  but  the  Dp/Hp  ^ata 
deviate  considerably  from  an  Arrhenius  plot  at  77  K,  at  which  temperature 
vheir  deduced  rate  coefficients  arc  about  the  same  as  at  173  K.  Furthermore, 
their  173  K result  differs  from  the  same  ratio  measured  by  Pcrsky  (1973)  by 
almost  a factor  of  2.  Only  the  room-temperature  measurement  appears  to  be 
reliable,  li  this  result  is  combined  with  their  Dp/CH^  expression,  a value 
for  thr  ratio  of  kp+j-ip/kp+d^  is  °^taincd*  In  Fig.  both  their  deduced 
rate  coefficient  and  the  result  of  ignoring  the  low-tcmperaturc  data  are  given. 


Foon  and  Reid  (1971)  conducted  an  extensive  scries  of  competitive 
fluorination  measurements  with  Hp  and  a variety  of  hydroca  rbons . They  used 
a packed,  heated  vessel  to  measure  the  consumption  of  reactants  bv  gas  chro- 
matography. They  found  that  no  photo^initiation  was  required,  the  F atoms 
being  produced  in  sufficient  quantities  by  the  equilibrium  dissociation  of  Fp. 
They  measured  the  rate  of  reaction  with  Hp  relative  to  CH^  and  CpH^  in  sepa- 
rate experiments.  A third  measurement  of  the  relative  rates  with  CH^  and 
C^ll  provided  an  additional  check  on  the  internal  consistency  of  the  measure- 
ments.  However,  in  spite  of  the  internal  agreement  demonstrated,  it  is  not 
obvious  how  the  complication  of  chain  reaction  (H  + Fp)  is  avoided. 

Jonathan  et  al.  (i°7icl  produced  F atoms  by  a discharge  through  CF^, 
then  observed  the  infrared  chemiluminescence  from  vibrationally  excited  HF 


k , cc/mol  sec 


1DOQ 


500 

X 


T.  K 

333  ?50 

—i r~ 


200 

X 


10 


14  _ 


10 


13  _ 


10 


12  _ 


MAYER  et  al  ‘BEBO  calculation]  119661 

KAPRALOVA.  MARGOLIN.  AND  CHAIKIN  (19701 

IGOSHIN  et  al.  1 19741 

HOMANN  etal.  119701 

fOON  AND  REID  119711 

KCMPAANO  WANNER  1 19721 
o OODONOV  st all  19711 

o ClYNE  stall  19731 

' HON  etal.119731 

WILLIAMS  ANO  ROWLANO  119731 
' RABIOEAU,  HECHT,  ANO  LEWIS  1 19721 


10 


11 


3 4 5 

1Q00/T,  1/K 


Fig.  4.  Rate  Coefficient  k~  for  F + H2‘ 
HF  + H Versus  lO^/T 


26 


4) 

u 

c 

it 

V- 

<y 

w. 

<y 

cd 


tT* 

* i 


o 

& 


41 

P 

.? 

e 

u 

CJ 

H 


c 

£ 

SJ 

U 

a 

rt 

o 

2 

> 


c£ 


T? 

u 

TJ 

h 

C 

■rf"‘ 

C 

rt 

X 

u 

h 

l* 

r- 

r- 

O" 

03 

pH 

S 

O 

c: 

X 

2- 

“0 

cc: 

TJ 

T3  _ 

■0 

c 

nJ 

m C 

flj  \S 
> :- 

V 

c 

C m 
« ^ 

u 

0 « 

C 

a 

P 

i — 
O 0s 
»-.  ^ 
4.  P' 

2r 

rl  x 

2 0 
s-o 

x s 

c: 

0 

0 

U, 

pC 

4J 

rt 

c 

O 

h 

1-1 

.2- 

s «- 
Zi  O' 

x 


aT  g 

U rn 


u 


CJ 

rS 


o 


>* 

~Z  a 
JT  o 
— a 

U » 5 

. * wv 

- rt  _ 

g 

X ^ CL 
. - Xi  «« 

r- 1 
<1 

£ 

u L. 

o o 
Cl  u 

f£*  2 

Ui  ^ 


U 

o 


*2 


u j? 
v 2 

^ x* 

*u  . — 
+-* 

Q — 

7:  .£ 


* i 

u « « 

“ » ? 

Ti  E K 


v 2 

u *_, 

15.5 

P“* 

n 


CJ 

u 

w 1 tn 

• F 

^ 2 X — 

—v  O 

C C7  w u* 

p W 

JL  cj 
X c 

®os°- 
2 b ii  ^ 0 

3^  “«i 

u ,b 

*3  y J O N 

*-  i 

X- 

2 « = 

V 
bO 

b *£ 

3 M> 

J*-  3 

V o * 

cfu 


NT 

v4 

nj 

TT 

’T 

• 

% 

% 

■ 

O 

0 

O 

0 

H 

Ui 

X 

p£ 

H 

a: 

0 

0 

'*ss. 

0 

00 

0 

cf\ 

0 

rn 

1 

vO 

t 

"0. 



J, 

X 

a 

(X 

V 

X 

X 

in 

4J 

r- 

<M 

U 
d 
o 
■— < 
u 

7: 


*r 

CO 


0 

O 

c 

0 

X 

CJ 

oj  no 

p-* 

u 

<m  > jr 

- > X 

rt 

X 

OCX 

* vO 

<M  1m  (NJ 

c? 

ci 

w 

« O 

*c*  ^ 
X 

X 

X jrrt  0 

O 

rs 

c £ 2 

u u 

c 

vu  u 2 

•* 

m 

r"' 

OC  m 

nj 

■^S 

tt  rf 

*r 

tN) 

m m 

0 

0 

» 

1 

1 1 

0 

0 

* 

JC 

m 

m m 

m 

ro 

H 

r- 

*p  r- 

? 

ni 

— 

r\l  N 

Z7 


Ignores  data  at  77  K. 


I 


produced  by  reaction  with  various  reagents,  including  H2  and  CH^ . From  the 
relative  intensities,  they  deduced  relative  F-atom  rate  coefficients. 

1 8 

Williams  and  Rowland  (197  3)  formed  F by  nuclear  reaction;  the  trans- 
lationaily  hoi  atoms  were  moderated  by  collisional  deactivation,  then  allowed 
to  react  with  pairs  of  reactants,  one  of  which  was  always  Products  were 

analyzed  by  iadio  gas  chromatography.  However,  concentrations  were  never 
measured  directly,  but  only  inferred  from  the  diminution  in  CH->CHF  yield. 

Results  of  the  above  competitive  studies  arc  listed  in  Table  3.  From 

Tables  2 and  4 can  be  calculated  k„Ltf  /k_  -tI  from  those  studies  in  which 

F+H  2 F+GM  4 

the  two  reactions  were  studied  separately.  With  the  exception  of  the  result 

of  Kapralova,  Margolin,  and  Chaikin  (1970),  the  various  values  for 

kr,rn  fall  between  0.  34  and  0.  88  at  room  temperature,  A possible  expla- 
r 

nation  for  the  high  value  obtained  by  Kompa  and  Wanner  (1972)  has  already 
been  discussed.  The  best  value  for  the  ratio  is  probably  0.  6 ± 0,  2. 

Measurements  of  the  F + CH^  reaction  rate  coefficient  arc  given  in 
Tabic  4.  The  absolute  studies  of  Kompa  and  Wanner  (1972);  Hon.  Axworthy. 
and  Schneider  ( 1973);  and  Clync,  McKenney,  and  Walker  (1973)  have  already 
been  discussed,  as  has  the  relative  measurement  of  Foon  and  Reid  (1971). 

Foon  and  Reid  studied  relative  fluorination  rates  for  several  pairs  of  hydro- 
carbons, ultimately  relying  on  the  assumed  rate  coefficient  for  the  F + C^Kjq 
reaction  arrived  at  by  Fettis  and  Knox  (1969)  after  a careful  review  of  all  the 
data  available  at  that  time.  Foon  and  Reid’s  (1971)  CH^  rate  coefficient  is 
about  60%  smaller  than  the  value  recommended  by  Fettis  and  Knox  (1969). 
However,  both  of  these  values  depend  on  the  assumption  that  the  activation 
energy  for  F + higher  hydrocarbons  is  0,  whereas  the  number  actually  may 
be  slightly  negative.  If  it  is  negative  by  a few  hundred  calorics,  then  the 
calculated  value  for  + would  be  larger  than  is  generally  quoted,  placing 
the  results  of  Fettis  and  Knox  (1969)  and  of  Foon  and  Reid  (1971)  in  better 
agreement  with  the  other  values  in  Table  4. 

One  other  absolute  measurement  of  F + CH^  has  been  reported  by 
Wagner,  Warnatz,  and  Zctzsch  (1971).  They  produced  F atoms  by  the  reaction 
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of  N f NF?  in  the  same  apparatus  used  by  Homann  ct  al.  (1970)  for  the  study 
of  the  F + reaction.  Their  result  of  4 . 9 X 10  cc/mol-sec  is  about  15% 
larger  than  that  of  Kompa  and  Wanner  (1972). 

The  measurement  of  Pollack  and  Jones  ( 1973),  which  is  somewhat  larger 

than  any  of  the  others  (Table  5),  was  made  relative  to  the  rate  of  the  reaction 

of  F + NO  + M = FNO¥  + M,  a reaction  that  docs  not  seem  to  be  sufficiently 

well  understood  at  present  to  serve  as  a calibration  standard  for  kinetics 

measurements.  Thus,  the  best  value  for  the  F + CH4  reaction  rate  co- 

13  * 

efficient  at  room  temperature  is  probably  4 t 1 X 10  cc/mol-sec,  making 
the  best  estimate  of  k^  2.6  - 1.4  x 10  slightly  larger  than  the  results  of 
the  direct  study. 

c . Complex  Systems 

Almost  a decade  ago,  Levy  and  Copclcnd  (1965)  attempted  to  deduce  the 
value  of  k^  from  the  inhibiting  effect  of  oxygen  on  the  reaction  of  with  F 2* 
They  believed  the  inhibition  was  due  to  one  of  two  reactions,  either 

H +Oz  + M -**  HO^  +M 

in  competition  with  the  chain  propagation  step 

H + F2—  HF  + F 


or  by  the  analogous 


F + O £ + M FOz  + M 

competing  with  Reaction  (4).  Assuming  the  oxygen  reacted  with  H atoms,  they 
deduced  an  activation  energy  for  Reaction  (4)  of  5 to  7 kcal/mol.  On  the  other 
hand,  the  assumption  that  inhibition  occurred  by  the  reaction  with  F atoms  im- 
plied an  activation  energy  for  Reaction  (5)  of  5 to  7 kcal/mol.  They  concluded,  on 
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the  basis  of  comparison  with  analogous  reactions  with  known  activation 
energies,  that  5 to  7 keal/mol  was  too  large  for  Reaction  (5),  but  reactions  of 
other  similar  atoms  and  radicals,  i»e»,  H,  D,  and  CH,  with  had  activa- 
tion energies  in  the  range  of  5 to  10  kcal/mol.  They  concluded,  therefore, 

chat  inhibition  was  by  reaction  with  F atoms,  and  Reaction  (4)  had  an  activa- 

1 2 

tion  energy  of  7 kcal/mol  with  a pre-exponential  factor  of  10  cc/mol-scc 
assumed.  These  calculations  were  made  by  taking  the  reaction  rate  coeffi- 
cient for  H + + M to  be  10^  * cc^/mol-sec,  and  that  for  F + + M to  be 

the  same.  However,  since  they  performed  their  work,  more  recent  experi- 
ments indicate  that  the  rate  coefficient  for  H + + M is  about  50  times 

larger  than  was  accepted  at  the  time  of  Levy  and  Copeland's  (1965)  work. 
Therefore,  if  the  competition  is  for  H or  F atoms,  the  activation  energy  for 
Reaction  (5)  or  (4)  would  be  1,8  to  3,7  kcal/mol.  Since  this  is  reasonable 
for  the  H + F^  reaction,  there  is  no  need  to  postulate  the  otherwise  unknown 
reaction  of  F + -f  M*  Therefore , Levy  and  Copeland's  (1965)  work  should 
not  be  taken  as  suggesting  anything  about  the  activation  energy  of  Reaction  (4). 
Later  (Levy  and  Copeland,  1968),  they  rejected  this  conclusion  on  other  grounds 

Rabideau,  Hecht,  and  Lewis  (1972)  used  a microwave  discharge  to  dis- 
sociate in  a teflon-coated  flow  system,  obtaining  about  66%  dissociation. 
The  H + H^>  in  a He  carrier,  was  mixed  with  a stream  of  F^,  And  the  reac- 
tion was  monitored  by  measuring  the  H and  F atoms  by  esr  spectroscopy. 
Since  it  was  recognized  that  both  reactions,  F + and  H + F^,  were 
occurring,  a computer  program  was  used  to  simulate  the  observed  atom  pro- 
files and  thereby  deduce  reaction  rate  coefficients.  At  room  temperature, 

12 

best  results  were  obtained  by  the  assumption  that  k . =■  4 0 1 X 10  and 
1 2 ** 

k^  = 2.  5 ± 0.  2 X 10  cc/mol-sec.  This  result  for  agrees  well  with  those 
of  foon  and  Reid  (1971)  and  of  Williams  and  Rowland  (1973)  but  is  a factor  of 
2.5  to  10  lower  than  most  of  the  other  determinations.  Rabideau,  Hecht,  and 
Lewis  (1972)  gave  no  indication  of  the  sensitivity  of  their  result  for  to 
changes  in  their  result  for  k^;  obviously,  a decrease  in  one  of  these  rate 
coefficients  would  induce  an  increase  in  the  calculated  value  for  the  other. 
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Igoshin,  Kulakov,  and  Nikitin  (1974)  determined  k^  from  the  stimulated 

emission  of  HF  molecules  formed  by  the  electric-discharge-initiated  reaction 

between  and  NF^  in  a laser  cavity  at  195  and  293  K.  A computer  model, 

which  took  into  account  initiation,  vibrational  energy  transfer,  spontaneous 

and  stimulated  emission,  and  wall  deactivation,  was  used  to  match  the  power 

output  from  the  lasing  medium.  The  accuracy  of  this  technique  must  depend 

on  the  values  assumed  for  the  other  reaction  rate  coefficients;  however, 

Igoshin,  Kulakov,  and  Nikitin  (1974)  do  not  give  detailed  values  for  individual 

rate  processes,  making  it  difficult  to  assess  the  validity  of  the  calculations. 

From  their  two  temperature  determinations,  they  calculated  an  Arrhenius 

expression  of  k.  = 9.3  x 10  ^ exp[  - ( 1 080  - 170)/RT]  cc  /mol  - sec.  Their 
* 13 

room- temperature  rate  is  2 X 10  cc/mol-sec,  almost  a factor  of  2 smaller 
than  that  of  Kompa  and  Wanner  (1972),  the  largest  value  reported.  Igoshin, 
Kulakov,  and  Nikitin  (1974)  separately  determined  the  rate  coefficient  for  the 
F + reaction,  which  provides  an  independent  check  on  the  reliability  of 
their  method.  They  obtained,  for  + 4.9  X 10  cxp[-(790  - 180)/RT] 

cc/mol-sec,  giving  k^^/k^  = * * exp(-290/RT).  Two  recent  direct  experi- 

mental measurements  of  this  ratio  jive  values  very  close  to  1.06  exp(360/RT) 
or  1.94  at  300  K.  The  value  deduce  f by  Igoshin,  Kulakov,  and  Nikitin  (1974) 
disagrees  with  these  more  reliable  results  by  approximately  a factor  of  2. 

d.  Theoretical  Calculations 

Wilkins  (1971  and  1972);  Jaffc  and  Anderson  (1971);  Jaffe,  Henry,  and 
Anderson  (1973);  and  Muckerman  (1971,  1972a,  1972b)  have  reported  the 
results  of  high-speed,  computer-generated  trajectory  calculations  of  the  over- 
all rate  coefficient  for  Reaction  (4);  however,  these  results  depend  sensitively 
on  the  parameters  of  the  potential  energy  surface  used,  and  since  there  are 
as  yet  no  precise  a priori  methods  for  generating  such  surfaces,  the  technique 
is  of  little  value  in  generating  independent  Arrhenius  factors  for  k^.  For 
determining  rate  coefficients  for  the  production  of  individual  vibrational- 
rotational  energy  states  of  product  molecules  or  for  examining  the  effect  of 
excess  reagent  internal  energy,  however,  such  techniques  do  provide  rough 
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information  that  might  not  be  available  otherwise.  Such  work  is  discussed 
in  more  detail  later  in  this  review.  Before  the  machine  calculations  were 
made,  several  estimates  and  calculations  of  k^  were  published.  Only  the 
BEBO  calculations  of  Mayer,  Schieler,  and  Johnston  (1966)  and  Mayer  and 
Schieler  (1968)  are  of  interest  here.  They  predicted  a room-temperature 
value  for  k^  of  about  6 X 10  , over  an  order  of  magnitude  smaller  than  the 

average  experimental  value. 


Two  other  theoretical  studies  are  of  some  interest.  Tully  (1974)  studied 

2 

the  role  of  spin-orbit  effects  on  Reaction  (4).  He  found  F(  Pl/2)  to  react 
almost  an  order  of  magnitude  slower  with  than  the  ground  state  Ff^P^^). 
Klein  and  Persky  (1 974)  studied  the  reaction  with  para -enriched  H^  and  con- 
cluded that  the  rotational  excitation  of  the  molecules  slews  the  reaction, 
but  the  effect  is  small,  and  important  only  at  very  low  temperatures. 

e.  Conclusions 

The  five  direct,  absolute  measurements  suggest  a room  temperature 
value  for  k^  of  1.5  - 0.5  x 10 1 ^ cc/mo l-sec.  On)y  one,  that  of  Homann  et  al. 
(1970),  was  obtained  for  a range  of  temperatures,  making  it  possible  to  cal- 
culate an  activation  energy  for  the  reaction.  The  room-temperature  value  of 

i 3 

Homann  et  al.  (1970)  is  1.0  x 10  ; the  competitive  studies,  together  with  an 

evaluation  of  the  F + CH^  rate  coefficient,  suggests  a slightly  larger  value 
of  k^,  but  with  correspondingly  larger  uncertainty.  A brief  analysis  of  the 
technique  of  Homann  et  al.  (1970)  suggests  only  possible  errors  that  would 
have  led  them  to  deduce  too  large  and  not  too  small  a value  of  k^. 

2.  VIBRATIONAL  PRODUCT  DISTRIBUTION 


Three  laboratories  have  studied  the  product  distribution  of  Reaction  (4). 
Pimentel  and  his  coworkers  [Kompa  and  Pimentel  (1968),  Parker  and 
Pimentel  (1969),  Krogh  and  Pimentel  (1972)]  examined  the  reaction  by 
measuring  gains  of  appropriate  spectral  lines  in  a chemical  laser.  Polanyi 
ar.d  coworkers  f Polanyi  and  Tardy  (1969),  Anlauf  et  al,  (1970),  Polanyi  and 
Woodall  (1972)]  and  Jonathan  et  al.  (1971c)  used  a discharge  flow  system  and 
made  measurements  in  a nonlasing  medium.  The  various  results  are  given  in 
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Table  5.  The  three  groups  have  converged  on  the  value  k^(v  = l)/k^(v  2) 

- 0.3.  The  two  groups  that  have  measured  k^(v  - 3)  and  k^(v  - 2)  now  sub- 
stantially agree  as  well  with  the  probable  value  for  the  ratio  of  0.48  to  0.5. 
There  is  one  recent  report  on  the  experimentally  measured  temperature 
dependence  of  the  relative  pumping  rates.  Coombe  and  Pimentel  (1973). 
using  laser  techniques,  found  k^(v  = 2)/k^(v  = 1)  r 2.  14  exp(254/RT)  and 
k4iv  - 3)/k4(v  = 2)  = 0.39  exp(117/RT). 

Because  HF  concentrations  are  determined  by  emission  measurements, 
the  formation  rate  of  HF(v  = 0)  has  not  been  ascertained  in  any  of  the  experi- 
ments. Rece  nt  Monte  Car  lo  computer  calculations  by  Wilkins  ( 197  2)(Fig.  5)  are  of 
interest.  His  results  suggest  that  k^(v  = 0)  = 0 and  that  the  activation  ener- 
gies for  v = 1,  2.  and  3,  are,  respectively,  1.5,  1.23,  and  1.25  keal/mol, 
implying  a difference  in  activation  energies  between  k^2)  and  k^(l)  of  270 
cal/rrrjl,  which  agrees  with  experiment.  Since  the  differences  in  activation 
energies  are  rather  slight,  we  prefer  at  this  time  to  assume,  for  computa- 
tional convenience,  the  same  temperature  dependence  for  each  of  the  k so 
that  each  k^(v)  can  be  expressed  as where  l_  is  temperature  in- 
dependent. W<t  accept  Wilkins'  (1972)  result  that  negligible  pumping  into  the 
zeroth  vibrational  level  takes  place. 

Recently,  Coombe  and  Pimentel  ( 1973b)  reported  the  effect  of 
rotational  energy  on  the  vibrational  energy  distribution  in  reaction.  They 
found  a greater  tendency  to  inversion,  especially  on  the  2 1 transition  when 

the  was  par  a -enriched  (85%  para). 

Anderson  (1970)  published  in  brief  the  results  of  a classical  trajectory 
study  of  the  reaction  between  H and  HF(v)  for  v = 0.  . . 5,  showing  the  threshold 
translational  energy  required  for  reaction.  Wilkins  ( 1973b),  with  a Monte 
Carlo  computer  technique,  calculated  the  rates  for  reactions  of  the  form 
I!  + HF(v)  F + H^v1)  for  v ^ 3,  4,  5,  6 and  v1  - 0.  1,  2.  His  results  suggest 


34 


k4(v) 


Table  5.  Relative  Pumping  Rates  for  F + — • HF(v)  + H 


k4(v)  (300  K) 

Technique 

Reference 

V = 0 

V = 1 

v = 2 

v = 3 

-- 

1.5 

. 1 

-- 

Laser  experiment 

Kompa  and  Pimentel 
(1968) 

-- 

0.  18 

1 

<1.  33 

Laser  experiment 
(T  = 539  K) 

Parker  and  Pimentel 
(1969) 

-- 

~0.  3 

1 

-- 

Laser  experiment 

Krogh  and  Pimentel 
(1972) 

-- 

<0.29 

1 

>0.47 

Discharge  flow 

Polanyi  and  Tardy 
(1969) 

-- 

0.  31 

1 

0.48 

Discharge  flow 

Anlauf  et  al.  (1970) 

0.31 

1 

0.47 

Discharge  flow 

Polanyi  and  Woodall 
(1972) 

- - 

0.29 

1 

0.76 

Discharge  flow 

Jonathan,  Melliar- 
Smith,  and  Slater 
(1971) 

-- 

0.  30 

1 

0.  5 

Discharge  flow 

Jonathan  et  al.  (1971c) 

-0 

0.  15 

1 

0.49 

♦ 

Semi-empirical 
calculation /Monte  Carlo 

Wilkins  (1972) 

0 

0.  53 

1 

0.  18 

Semi -classical 
calculation /Monte  Carlo 

Jaffe  and  Anderson 
(1972) 

0 

1. 28 

1 

0. 0023 

Variational  theory/ 
Monte  Carlo 

Jaffe,  Henry  and 
Anderson  (1973) 
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r*  - 


Fig.  5.  Values  of  k4(v)  Versus  103/T  for  F + H2  — 
HF(v)  + H for  Various  v as  Computed 
by  Wilkins 
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that  all  of  these  reactions  proceed  with  activation  energies  of  approximately 
500  to  600  cal/mol  and  that  at  high  temperatures  total  probability  for  reaction 
of  a given  v into  all  possible  v'  is  about  0.  1 .o  0.  2;  the  rate  coefficients  are 
tabulated  in  Table  6.  The  results  can  be  simplified  and  summarized  approx- 
imately (within  20  to  30%)  by  the  following  expressions  for  various  k 

k_4(6,v*)  = f6(v')  (1.9  X 1013)  exp(-  580/RT) 
f6(0)  = 0.22,  f6(l)  = 0.22,  f6(2)  = 0.56 

k_4(5,v')  = f5(v')  (1.  I X 1 0*  3)  cxp(- 5 1 0/RT) 
f5(0)  = 0.36,  f5(l)  = 0.64,  f 5 ( Z)  = 0 

k_4(4,v')  = f4(v’)  (7.4  X 1012)  exp(-460/RT) 
f4(0)  = 0.5,  f4(l)  = 0.5,  f4(2)  = 0 

k_4(3,v')  = f3(v')  (1.62  X 1013T"0,01)  exp(- 835 /RT) 

f3(0)  = 1 


B.  H 4-  F,  — HF  4 F 

In  comparison  with  the  previous  reaction,  experimental  work  on  the 
reaction 


H + F2-HF  + F (5) 

has  been  meager.  The  data  for  k^,  both  experimental  and  theoretical,  are 
tabulated  in  Table  7.  There  has  been  only  one  experimental  determination 

of  kj.  over  a temperature  range  sufficient,  to  permit  calculation  of  Arrhenius 

. , 14 

rate  parameters.  Albright  ct  al.  (19o9)  obtained  a value  of  k,.  =■  1.2  X 10 
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k *#  v* ) 

Table*  6.  Computed  Rate  Coefticients  for  H + HF(v)  <,  H^fv')  + Fa 


V 

v' 

k = ATne'E/RT 

loii  A 

n 

E,  cal/mol 

6 

0 

.2.66 

-0.  03 

570 

1 

12.41 

0.085 

540 

2 

13.02 

0.01 

565 

5 

0 

12.  £4 

0.015 

510 

1 

12.92 

-0.005 

575 

4 

0 

12.82 

-0.44 

570 

i 

12.  10 

0.  13 

380 

3 

0 

13.21 

-0.01 

835 

aFroin  Wilkins  (1973b) 
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exp(-2400/ R T)  in  a discharge  flow  system  at  temperatures  of  294  to  565  K, 
using  mass  spectrometry  to  monitor  species.  However,  their  determination 
of  the  rate  coefficient  for  H + C^,  obtained  in  the  same  apparatus  at  the 
same  time,  is  larger  than  some  oth*  r evaluations  by  an  order  of  magnitude 
at  300  K;  therefore,  their  value  of  is  somewhat  subject  to  •*  It  by 
association."  The  measurement  of  Rabidcau  (1972)  was  made  in  an 
system  in  which  both  Reactions  (4)  and  (5)  were  occurring.  Observed  F- 
and  H-atom  concentrations  were  matched  by  computer  calculations;  the  best 
fit  was  obtained  with  their  resultant  values  of  and  k,.,  which,  possibly, 
are  subject  to  greater  uncertainty  than  their  quoted  25%  (k^)  and  8%  (k^h 

One  earlier,  indirect  experimental  measurement  has  also  beer  reported. 

Levy  and  Copeland  (1968)  examined  the  kinetics  of  the  H2-F2  photochemically 

initiated  reaction  at  288  K.  In  the  presence  of  O^i  inhibition  was  noted,  which 

was  accounted  for  by  the  reaction  H + O^  + M,  where  in  their  system  M = N^i 

principally.  The  value  of  k^  was  determined  from  the  ratio  of  rates  of  reaction 

for  H + to  H + O^  + N7.  Assuming  a value  of  9.  2 x 10*^  cm^-moT  ^ sec  * 

for  the  latter  rate  coefficient,  Levy  and  Copeland  (1968)  calculated  the  value 

1 2 

of  k-  at  288  K to  be  1,8  x 10  . Then,  assuming  a pre-exponential  value  of 

9 ^1 

10  * , they  calculated  an  activation  energy  of  1.5  kcal/mol.  The  rate  coef- 

ficient for  the  competitive  H + C>2  + N2  reaction  is  not  accurately  known;  at 
higher  temperatures,  N2  seems  to  be  slightly  less  efficient  than  Ar  (Raulch 
(1972)],  suggesting  a value  of  approximately  5 x 10*^  cm^-mol  ^ sec  ^ 
for  the  three  body  process,  so  that  k^  would  be  approximately  0.  5 times  the 
value  t.nat  they  obtained  and  the  activation  energy  would  be  slightly  larger. 

Two  theoretical  determinations  are  also  listed  in  Table  7.  The  BF130 

calculation  by  Mayer,  Schieler,  and  Johnston  (1967)  predicts  a room- 

temperature  value  for  k^  that  is  almost  certainly  too  large  by  about  an  order 

of  magnitude . Wilkins  ( 197  3a)  obtained  a value  for  kj.  of  T ^ ^ 

14  04  ^ 

exp( -2840/RT),  or  approximately  10  * exp(- 2450/RT)  between  200  and 

500  K,  by  three-dimensional  classical  Monte  Ci*rlo  trajectory  calculations, 
yielding  a room-temperature  value  about  25%  smaller  than  that  of  Albright 
et  al.  (1969). 
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It  is  evident  that  another  determination  of  k^  is  needed;  such  an 
experiment  is  under  way  in  this  laboratory.  In  its  absence,  however,  the 
value  of  Albright  et  al.  (1969)  is  recommended,  with  the  cautionary  note 
that  it  may  be  somewhat  too  large. 

The  energy  distribution  of  Reaction  (5)  has  been  studied  in  two  lab- 
oratories. Jonathan,  Melliar-Sm ith , and  Slater  (1970,  19?!a,  1972)  have 
reported  successively  improved  measurements  of  the  pumping  distribution 
in  a discharge  flow  system.  Polanyi  and  Sloan  (1972),  using  a somewhat 
different  discharge  flow  technique,  have  also  published  distribution  numbers. 
Roth  studies  report  measurements  of  emission  from  v = 1 through  9;  Jonathan’s 
(1972)  work  includes  data  on  the  relative  rate  of  formation  of  the  v = 0 
level  monitored  by  absorption  techniques.  Jonathan  (197  2)  used  the  measured- 
relaxation  method,  while  Polanyi  and  Sloan  (1972)  used  the  a rrested- relaxation 
technique.  In  principle,  the  latter  should  be  the  more  direct  technique,  but, 
in  practice,  both  procedures  require  great  care.  In  both  laboratories,  the 
emission  intensities  were  converted  to  population  densities  by  means  of 
Einstein  coefficients  calculated  by  the  method  of  Cashion  (1963).  Herbelin 
and  Emanuel  (1974)  have  recalculated  the  Einstein  coefficients  and  have  ad- 
justed the  data  of  Jonathan,  Melliar-Smith,  and  Slater  (1970,  1971a,  1972) 
and  of  Polanyi  and  Sloan  (1973)  to  produce  revised  sets  of  rate  coefficients. 

The  original  and  the  corrected  results  are  shown  in  Table  8.  The  g(v)  shown 
are  the  relative  k's  normalized  to  sum  to  unity.  Thus,  k. ^ ( v ) = g(v)k^  (total). 

In  the  original  work,  the  k’s  were  normalized  by  arbitrarily  setting  kmax  = 1. 

Wilkins'  (1973a)  classical  Monte  Carlo  Trajectory  calculations  (Fig.  6) 
of  the  distribution  of  vibrationally  excited  HF  product  molecules  in  Reaction  (5) 
indicate  that  no  molecules  are  formed  in  the  levels  with  v £2  or  > 7.  The 
giv>  obtained  at  3t0  K from  his  calculations  are  g(3).  . .g(6)  0.  10:0.  23:0.43: 

0.  25,  suggesting  that  the  low  vibrational  levels  observed  in  the  experiments 
are  formed  by  vibrati onal- vibrational  (V.V)  processes  and  from  the  "cold*' 
Reaction  (4)  and  the  very  high  vibrational  levels  are  populated  by  VV  processes, 
which  are  expected  to  be  very  rapid  for  large  v.  It  is,  therefore,  assumed 
that  reaction  takes  place  into  levels  v = 4 through  8 only  and  Polanyi  and 
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T.K 

1000  000  333  250  200 


Fig.  6.  Values  of  k^(v)  Versus  1 /T  for  H + F^  -*HF(v)  + F 
for  Various  v as  Computed  by  Wilkins  (1973a) 
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Corrected  with  Einstein  coefficients  of  Herbelin  and  Emanuel  (1974)* 


of  Herbelin  and  Emanuel  (1974).  are  renormalized  so  that  g^(v)  - 1.  then 
the  values  g^(3):g5(4);g^(  5):g^(6):  = 0.07:0.14:0.35:0.44  arc  obtained.  These 
numbers  arc  expected  io  be  near  the  values  that  will  bo  agreed  upon  eventu- 
ally and  arc  recommended  for  computational  purposes. 


44 


IV.  ENERGY-TRANSFER  PROCESSES 


A.  VIBRATION  A E-TRANS1  J\TIONA  L (ROTATIONAL) 

ENERGY  TRANSFER  INVOLVING  HF 

Vibrational -translational  ( rotational)  t VTR)  energy-transfer  processes 
drive  the  vibrational  populations  toward  equilibrium  with  the  translational 
degrees  of  freedom  with  a resulting  loss  of  potential  power  in  a laser  system. 

The  relaxation  processes  for  the  hydrogen  halides  have  posed  a chailcng 
to  both  experimentalists  and  theoreticians,  Tnese  molecules  provide  a sensi- 
tive test  of  any  theory  of  vibrational  energy  transfer  because  of  their  large 
rotational  velocities,  their  strong  intermolccular  attractive  forces,  and  the 
large  effect  of  isotopic  substitution.  Theoretical  predictions  based  on  the 
Schwartz,  Slawsky,  and  Herzfeld  (1(?/  (SSH)  vibrational-translational  (VT) 
theory,  the  vibrational-rotationr.l  (VR)  theorv  of  Moore  ( 1963),  or  Millikan 
and  White's  ( 1963)  semi-empirical  correlation  do  not  reproduce  the  experi- 
mental data  for  the  hydrogen  halides.  A VR  model,  developed  by  Shin 
(1971b,  1973),  predicts  results  in  reasonable  agreement  with  the  hydrogen 
halide  self-relaxation  data.  In  this  model,  based  on  a collisior  between  a 
rotation-averaged  oscillator  and  a rigid  rotor,  the  two  H atoms  interact 
strongly  in  close-in  collisions  to  transfer  vibrational  and  rotational  energies. 
Classical  trajectory  calculations,  as  discussed  by  Wilkins  (1974),  have  been 
performed  in  only  a few  cases.  Because  of  the  difficulties  of  the  theoretical 
calculations,  only  experimental  data  exist  for  most  of  the  relaxation  proces- 
ses. For  some  of  the  processes,  there  are  neither  experimental  data  nor 
theoretical  calculations. 

The  VTR  deactivation  of  the  HF  can  he  described  by 

Mi 

^6(v,  v-Av) 

I!F(v)  + M.  1 HF( v - ^v)  + M. 
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M 

where  k,,1  A k is  the  rate  coefficient  of  the  reaction,  v represents  the  initial 
6(v,  v-Av)  r 

vibrational  level  of  HF,  and  Av  is  the  number  of  quanta  transferred  io  rotational- 
translational  energy.  Relaxation  rate  coefficients  have  been  measured  for  the 
self-deactivation  of  HF  and  for  its  deactivation  by  other  hydrogen  halides,  other 
diatomic  molecules,  and  some  atoms  and  polyatomic  molecules.  Most  of  the 
experimental  measurements  have  been  made  for  v = Av  = 1;  only  a few  studies 
have  reported  rates  for  the  higher  vibrational  levels. 

The  complexity  of  Shin's  (197  1b,  1973)  theory  and  the  specific  nature  of 
the  classical  trajectory  calculations  preclude  any  general  conclusions  on  the 
v and  Av  dependence  of  the  reaction  rates.  Since  the  dcactivations  of  the  upper 
vibrational  levels  affect  the  power  of  a laser,  kinetic  models  of  the  laser 
medium  must  include  these  upper  levels.  The  harmonic  oscillator  (HO)  model 
predicts  the  rate  coefficients  for  the  higher  vibrational  levels  to  bo  given  by 
k ( v v t)  = q)  for  Av  = 1 and  the  rr*"!tiquantum  transitions,  Av  > 1,  to  be 

negligible  compared  to  single  quantum  transitions.  HF  is,  however,  far  from 
an  ideal  harmonic  oscillator,  and  the  HO  model  cannot  be  expected  to  give  a 
precise  description  of  the  relaxation  processes.  The  extent  to  which  anharmon- 
icity  changes  the  v dependence  of  the  HO  predicted  rates  depends  on  the 
collision  model  assumed.  It  is  not  clear  whether  the  v dependence  of  k in- 
creases or  decreases  as  a result  of  anharmonicity . The  use  of  the  Landau- 
Teller  approximation  for  the  calculation  of  the  relaxation  rates  in  the  harmonic 
oscillator  model  results  in  a negligible  probability  for  multiquantutn  exchanges. 
This  approximation  is  based  on  a consideration  of  only  the  strong  repulsive 
intermolecular  forces  and  may  break  down  completely  where  .he  collision 
partners  exert  strong  attractive  forces.  Trajectory  calculations  have  been 
performed  for  a few  deactivation  processes,  and  these  results  are  used  as  a 
guide  for  estimating  the  higher  level  rate  coefficients  of  those  few  processes. 

In  tli i a review,  most  extrapolations  of  experimental  data  for  the  deactivation 
rates  to  higher  vibrational  levels  will  be  hased  on  the  IlO  model  for  lack  of  any 
other  information.  The  v dependence  as  well  as  the  multiquantum  transitions 
remain  the  largest  uncertainties  in  kinetic  models  of  chemical  lasers  and 
certainly  will  need  revision  when  critical  experiments  can  be  devised  to  measure 
them . 
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Most  experimental  investigations  result  in  the  measurement  of  a 
relaxation  time  t at  a given  pressure  p and  temperature  T.  To  convert  this 
relaxation  time  to  a rate  coefficient  k in  units  of  cc/mol-scc,  certain  assump- 
tions concerning  the  nature  of  the  relaxation  process  must  be  made.  According 
to  the  HO  model, 


micity  of  the  reaction,  and  R is  the  universal  gas  constant. 

In  the  following  paragraphs,  different  chaperones  arc  treated.  An 
analytical  expression  hag  been  fitted  to  the  results  of  each  experimental 
study  in  which  the  deactivation  reactions  were  studied  over  a range  of  tem- 
peratures. These  expressions  tit  the  data  with  an  accuracy  of  10  to  15%. 
Most  of  the  results  are  summarized  in  Tables  9 through  12. 


Self- relaxation  of  HF  has  been  studied  at  the  higher  temperatures  (above 
1 350  K)  in  shock  tubes  and  at  lower  temperatures  by  laser-induced  fluorescence 
studies.  The  laser-induced  fluorescence  technique  has  been  used  mostly  at 


reflected  shock  waves  have  yielded  data  over  a wide  range  of  temperatures. 

The  relaxation  of  HF  has  been  studied  more  extensively  than  any  of  the  other 
hydrogen  halides,  with  rates  reported  between  295  and  5000  K. 

Shock-tube  studies  of  Solomon  et  al.  (1971),  Bott  and  Cohen  (197  1a),  Just 
and  Rimpcl  ( 1972),  and  Vas  il'ev  et  al . gave  values  of  pr  for  HF-HF  re- 

laxation at  1400  K between  0.  19  and  0.25  [iftcc/atm,  or  k between  4.7  and 
6.2  x 101  1 cc/moI-:;ec.  At  4000  K,  Solomon  ot  al.  ( I 97  1 ) measured  somewhat 
faster  rates  and  Just  and  Rimpcl  (1972)  measured  somewhat  slower  rates  than 
those  found  by  Bott  and  Cohen  (1971a)  and  by  Vasil'ev  et  al.  ( 1973). 


UT 


RT 

r r : 


(1,0) 


- k 


(0.  i)  k(l,  0)[1  - exp(-  AE/RT)] 


j 

i 


1 . M = HF 


room  temperature,  although  some  studies  in  heated  or  cooled  cells  or  behind 


47 


Table  9.  Experimental  Rates  for  VTR  Relaxation  of 
HF(v  - 1)  by  Various  Collision  Partners 


HF( 

V r 1 ) + M 

k.(Pr)-1 

HF(v  = 0)  + M 

M 

T,  K 

log  A 

E. 

cal/mol 

fV*; 

Reference 

Ar 

1330- 1580 

12.85 

29.800 

O 

o 

Johnson  and  liritlon  (1964)* 

Ar 

900-1900 

10.24 

14,000 

3.  6.  15 

Cglukian  ( 1964  )b 

Ar 

1300-1700 

12.49 

27,  300 

5 

Seery  and  Britton  (1966)C 

Ne 

1650-2700 

12.18 

23,900 

0.5,  1 

D.esen  (1966) 

Ne 

1400-2000 

13.30 

35,000 

0.5 

Diesen  (1968) 

Ar 

1200-1500 

13.57 

31,700 

10 

Just  and  Rimpel  (1970)^ 

Ar 

1400-2600 

13.55 

34,700 

10.  20 

Breshears  and  Bird  (1973) 

F2 

1400-2600 

13.99 

34.800 

10.  20 

Breshears  and  Bird  (1973) 

Ar 

1575-2540 

13. 10 

32,700 

0.5 

Mastromonico  (1974) 

A r 

-- 

12.66 

28.500 

(review  of 
da  a) 

Lloyd  (1971) 

* The  tabulated  results  are  for  the  six  runs  with  5%  K^.  1*  «•  irrVi  of  nine  runs  in 

10%  in  Ar  mixtures,  Johnson  and  Britton  (1964)  obtained  lot*  A 3 11.54  and 
E 3 20,680  cal/mol;  in  a s cries  of  five  runs  in  20%  in  Ar , they  obtained 
log  A 3 9.61  and  E : 11, 025  eal/mol.  They  considered  their  5%  results  to 
be  the  most  reliable. 

^The  results  are  based  primarily  on  the  twelve  6%  runs.  The  3%  runs  yielded  a 
considerably  smaller  value  for  E.  The  scatter  in  the  15%  runs  was  loo  large  to 
permit  a reliable  assignment  of  E. 

cResults  are  based  on  data  of  Johnson  and  Britton  (1964)  and  Seery  ar.d  Britton  (1966). 
In  mixtures  of  3%  F^,  20%  Hr,  and  75%  Ar,  Seery  and  Britton  obtained  log  A * 13,  15 
rnd  E = 21,  100:  in  10:20:70  mixtures,  they  obtained  log  A - 11.67  and  F s 19,600 
cal/mol. 

f*Just  and  Rimpel  ( 1 9 70 ) also  expressed  results  in  the  form  A^expf-D/RT), 

obtaining  log  As  11. 54,  ns -1.9,  and  D * 36,  700  cal/mol. 


43 


Table  10.  Calculated  Vibrational  Relaxation  Rate  Coefficients  of  HF(v)  by 
F Atoms 

HF(v)  + F 6(v'  yl)  HF(v')  + Fa 


kv  v,  - ATn  cxp( - E/RT) 
log  A n E,  cal/mol 


1 

0 

14 . 8_> 

-0.766 

3100 

2 

1,0 

14.79 

-0.708 

3075 

3 

2 

14.99 

-0.724 

364  5 

3 

1 

14.80 

-0.757 

3720 

6 

5 

15.71 

-0.712 

3645 

6 

4 

15.  56 

-0.741 

3640 

6 

3 

15.15 

-0.70 

3 59C 

6 

2 

14.89 

-0.658 

3810 

aFrom  Wilkins  (1973c).  Calculated  values  for  k.  , and  k,  , 

3,  v*  6f  v‘ 

include  comparable  contributions  from  both  reactive,  i.e., 
H-atom  transfer,  and  nonreactive  processes. 
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Tabic  11.  Experimental  Deactivation  Rate  Coefficients  of  HF(v)  by  H Atoms 


HF (v)  + H HF(v')  + H.  v’  < v 


V 

$ 

V 

to"11  k( 300  K). 
cc/mol-scc 

Reference 

1 

0 

1,4  ± 0.4 

Heidner  and  Bott 

1 

0 

<0.  1 

Quigley  and  Wolga  (1974) 

1 

0 

19 

Kwok  and  Wilkins  (1974) 

2 

l.o 

110 

Kwok  and  Wilkins  (1974) 

3 

2.  i.O 

180 

Kwok  and  Wilkins  (1974) 

1 

0 

4 

Kwok  and  Wilkins*5'  c 

2 

l.o 

40 

Kwok  and  Wilkins*5,  C 

3 

2.  i,0 

80 

Kwok  and  Wilkins*5*  C 

aR,  F.  Heidner  and  J.  F.  Bott,  J.  Chem.  Phys.  63, 

1810  (1975). 

^Recalculation  of  data  based  on  assumption  that  HF-HF  VT 
deactivation  scales  as  v^. 

M.  A.  Kwok  and  R.  L.  Wilkins,  The  Aerospace  Corporation, 
private  communication. 
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Tabic  12.  Computed  Rate  Coefficients  for  Deactivation  of  HF(v)  by  H Atoms1 

HF(v)  + H — (V'  V',)  HF(v')  + H 


V 

v1 

k^1  , = ATn 

6,  v,  v' 

cxp(-E/RT) 

log  A,  cc/mol-soc 

r* 

E. 

cal  /mol 

1 

0 

13.96 

-0.28 

1170 

2 

0 

13.63 

-0.08 

900 

2 

1 

12.64 

0.20 

7 20 

3 

0 

13.59 

-0.06 

760 

3 

1 

13.45 

-0.08 

680 

3 

2 

13.  57 

-0.09 

890 

4 

0 

14.  17 

-0.  32 

810 

4 

1 

12.82 

0.09 

640 

4 

2 

13. 10 

0.03 

670 

4 

3 

13.21 

0.0 

520 

5 

0 

13.98 

-0.23 

640 

5 

1 

13.60 

-0.  15 

610 

5 

2 

13.  30 

-0.04 

710 

5 

3 

13.50 

-0.06 

670 

5 

4 

13.34 

-0.05 

530 

6 

0 

13.88 

-0.  19 

620 

6 

1 

13.60 

0.  16 

570 

6 

2 

13.61 

-0.  15 

7 50 

6 

3 

13.23 

-0.01 

670 

6 

4 

13.82 

-0.21 

600 

6 

5 

13.96 

-0.22 

630 

From  Wilkins  (197  3b).  Values  for  n have  been  rounded  off 
to  two  significant  figures;  those  for  E to  the  nearest  10  cal. 
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The  laser -induced  fluorescence  technique  has  been  used  to  measure  HF 
relaxation  times  between  295  and  2400  K.  We  have  fitted  the  results  of  these 
various  determinations  with  function*  describing  the  temperature  dependencies 
and  listed  them  in  Table  9.  They  are  also  plotted  in  Fig.  (7).  The  results  of 
Lucht  and  Cool  (1974),  Stephens  and  Cool  (1972),  Hancock  and  Green  (1972a), 
Airey  and  Fried  (1971),  and  Ahl  and  Cool  (1973)  have  been  adjusted  for  the 
gasdynamic  correction  discussed  by  Bott  (1974a).  The  measured  relaxation 
times  agree  within  a factor  of  1 . 6 at  any  temperature  between  295  and  1000  K, 
with  the  exception  of  the  early  measurement  of  Airey  and  Fried  (1971)  at 

350  K,  which  is  somewhat  faster.  The  room -temperature  value  for  k is 

12 

1 X 10  cc/mol-sec,  within  about  20%.  The  high-temperatu ~e  data  combined 
with  those  at  low  temperature  show  that  the  HF-HF  VTR  relaxation  times  must 
have  a maximum  near  1200  K (Lucht  and  Cool  (1974),  Bott  (1972),  Osgood, 
Javan,  and  Sackett  (1972),  Fried,  Wilson,  and  Taylor  (1973),  Hinchen  (1973a, 
1973b),  and  Bott  and  Cohen  (1973b))  . This  maximum  wac  confirmed  by  Blair, 
Brcshears,  Schott  (1973)  in  an  experimental  study  covering  the  range  between 
600  and  2400  K.  However,  their  numerical  results  for  k,  in  the  neighborhood 
of  1500  K,  are  too  large  by  a factor  of  approximately  1.5.  The  phenomenon 
of  a relaxation  time  maximum  has  been  observed  or  predicted  for  other  dia- 
tomic molecules  (Donovan,  Husain,  and  Stevenson  (1970b)]  but  at  much  *ower 
temperatures.  The  increase  in  deactivation  rate  at  low  temperatures,  of 
course,  is  due  to  the  effects  of  an  intermolecular  attractive  force,  but  the 
origin  of  that  force  is  not  understood. 

Shin's  (1970,  1971b)  model  predicted  relaxation  rates  for  HC1,  HBr,  and 
HI  that  substantially  agree  with  experimental  data  above  700  K.  Applying  the  same 
theory  to  HF-HF  and  DF-DF  VR  energy  transfer.  Shin  (1971a,  1973)  found  it 
necessary  to  include  strong  dipoJe-dipole  attraction  to  reproduce  the  experi- 
mental data  between  300  and  4000  K. 

•c 

The  trajectory  calculations  by  Berend  and  Thommart  n for  HF-HF  and 
DF-DF  deactivation  gave  good  agreement  with  the  experimental  data  when  a 

G.  C.  Berend  and  R.  L.  Thommarson,  McDonnell-Douglas  Corporation, 
private  communication. 
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VASIL' EV  el  all  1973 1 


1 x 1013 


1 x 10 


12 


BOTT  AND  COHEN  (1971al 

JUST  AND  RIMPEL  (19721 
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HINCHEN  (19731 
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Experimental  Values  of  Rate  Coefficient  k^j 
HF(1)+  HF-HF(O)  + HF  (VT  Deactivation)  Versus  T 
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sufficient  number  of  J levels  were  considered.  Berend  and  Thommarson's 
(197  3a)  initial  calculations  were  based  on  fewer  J levels  than  were  necessary  to 
describe  adequately  the  relaxation  at  low  temperatures.  They  observed  that 
belosv  1000  K the  deactivation  resulted  from  the  collision  of  HF(v  = J)  with 
HF(v  - 0)  and  the  formation  of  internally  excited  (HF),  dimers.  The  excited 
dimer  lives  long  enough  to  allow  a redistribution  of  energy  so  that  neither  of 
the  HF  molecules  is  vibrationaliy  excited  when  the  dimer  dissociates. 

Both  HF  and  DF  self-relaxation  rates  have  been  described  by  two-term 

expressions  by  Bott  and  Cohen  ( 1973b)  that  fit  the  experimental  data  over  the 

temperature  range  of  295  to  4000  K.  At  temperatures  below  about  1000  K,  the 

low-temperature  term  dominates  the  high-temperature  term  and  Pt  is 

proportional  to  T^.  The  probability  p of  HF  self-deactivation  is  0.01  at  300  K 

and  extrapolates  to  about  0.02  at  200  K.  Low-temperature  relaxation  rates 

for  HF  as  yet  have  not  been  published;  however,  studies  by  Lucht  and  Cool 
- 1 3 

suggest  a T ‘ temperature  dependence  for  p below  300  K,  which  means 
- 0 8 < rt 

k « T * and  pT  oc  T . The  experimental  data  from  295  to  4000  K can  be 

fitted  by  the  expression  k = 3 x 10^T  * + 3.5  X 10^T^*^  cc/mol-sec,  where 

the  high-temperature  term  has  been  chosen  for  computational  ease  to  be  of  the 

form  T*1  instead  of  the  usual  Landau- Teller  exp(B/T)*^  dependence.  However, 

-0  8 

to  approximate  more  accurately  the  expected  T~  * dependence  at  low  tempera- 
tures, the  expression  k = 1 x i0**T  ® ^ should  be  used.  This  expression  is 
valid  for  temperature*  between  about  1000  and  200  K, 

The  question  of  how  to  scale  the  VT  self- relaxation  rate  with  v remains.  In 
past  reviews,  the  harmonic  oscillator  rule  v j)  ~ v^(i  0)  ^as  ^ecn  recom- 
mended  [Cohen  (1971,  1^72,  1974)].  However,  recent  efforts  to  fit  pulscd-laser 
data  as  well  as  flow-tube-kinetics  experiments  suggest  a strong  v dependence, 
with  k^ , ^ ^ ^K(t  0)*  uncertainty  in  the  higher  k^  ^ 1)  *ar£c*  kut  there 


R.  A.  Lucht  and  T.  A.  Cool,  Cornell  University,  private  communication. 

M.  Kwok,  The  Aerospace  Corporation,  private  communication,  and 

N.  Cohen,  The  Aerospace  Corporation,  unpublished  work. 
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is  some  indication  that  the  k's  increase  to  a maximum  around  v = 4.  There  is 
no  direct  evidence  for  multiquantum  transitions,  partly  because,  to  a good 
approximation,  any  multiquantum  transition  scheme  can  bo  closely  matched 
by  a single- quantum  scheme.  For  example,  allowing  all  multiquantum  trans- 
itions and  setting  = 1 0)  almost  the  same  deactivation  rate  for 

a given  v level  as  the  harmonic  oscillator  rule,  even  though  the  total  vibra- 
tional energy  of  the  system  would  be  removed  at  a faster  rate.  These  problems 
are  discussed  further  in  Section  111 .15. 

2.  M = Ar,  He 

Butt  and  Cohen  (1971a)  and  VasiFev  ct  at.  (1973)  have  reported  relaxation 
times  for  HF-Ar  relaxation  between  1500  and  5000  K that  substantially  agree 
with  each  other.  Blair,  Broshears,  and  Schott  (1973)  measured  somewhat  faster 
HF-Ar  relaxation  times  between  800  and  2400  K than  the  overlapping  and  extra- 
polated data  of  Bott  and  Cohen  (197  la)  and  VasiFev  et  al.  (1973).  Their  data 
also  show  some  curvature  on  the  Landau-Teller  plot,  whereas  those  of  Bott  and 
Cohen  (1971a)  and  Vasil'ev  et  al.  (1973)  do  not.  The  curvature  could  be  real  or 
the  result  of  the  effects  of  impurities  such  46  H^O  or  residual  HF,  which  become 
more  important  at  low  temperatures.  Bott  and  Cohen  (1971a)  found  He  to  be 
only  a factor  of  two  more  efficient  than  Ar  for  deactivating  H7\  Room-tempera- 
ture measurements  by  Airey  and  Fried  (1971),  Hancock  and  Green  (1972b),  and 
Hinchen  (1973a)  led  on’y  to  upper  limits  on  the  relaxation  rate  coefficient  of 
300  sec  * Torr  * at  350  K [Airey  aud  Fried  (1971)]  and  60  nec  * Torr  at  295  K 
[Hancock  and  Green  (1972b)  and  Hmchcn  (1973a)].  These  results  are  consistent 
with  extrapolations  from  the  shock-tube  data. 

The  SSH  VT  calculations  do  uot  adequately  predict  the  data  for  the  rare- 
gas  deactivations  of  the  hydrogen  halides  of  Bott  and  Cohen  (1973)  and  Zittel 
and  Moore  (1973),  The  trajectory  calculations  of  Bcrend  and  Tommarson 
( 1973b)  for  the  relaxation  of  HF  by  Ar  agree  within  a factor  of  three  with  the 
shotk-tube  data  but  do  not  have  as  steep  a temperature  dependence.  They 
concluded  that  a reduction  of  the  range  parameter  or^  ^ by  2u  to  30%  would 
lead  to  much  better  agreement  with  the  shock-tube  data.  Such  a variation 


is  well  within  the  range  of  uncertainty  of  a.  Variation  of  20  to  30%  in  the 
range  parameter  does  not  substantially  improve  the  SSH  predictions.  Steele 
and  Moore  ( 1974)  made  approximate  calculations  for  the  relaxation  of  HC1 
(DC1)  by  Ar  and  He  and  concluded  from  the  qualitative  fit  to  their  data  that 
both  rotational  and  translational  degress  of  freedom  play  important  roles. 
Experimental  data  and  theoretical  calculations  for  HF-Ar  vibrational  relaxa- 
tion are  shown  in  F;g.  8.  HF -He  data  [Bott  and  Cohen  (1971a),  Smith  (1974)] 
are  somewhat  more  sparse;  they  are  plotted  in  Fig.  8 also.  He  is  about  2 to  3 
times  more  efficient  than  Ar  throughout  the  temperature  range  of  the  measure- 
ments . 


3.  M = F 

Bott  and  Cohen  (1971b)  and  Solomon  et  al.  (1971)  observed  the  formation 
of  HF(v  = 1)  in  the  presence  of  F atoms  produced  by  dissociating  F^  in  one  case 
[Solomon  et  al.  (1971)]  and  SF^  in  the  other  case  [Bott  and  Cohen  (1971b)]  at 
temperatures  between  1900  and  3500  K in  shock  tubes.  Solomon  et  al.  (1971) 


reported  rate  coefficients  for  kf,.  of  5 X 10*^  cc/mol-scc  and  Bott  and 

h( 1 , 0) 

Cohen  (1971b),  3x10  ° cc /mol - sec  at  2500  K.  Blaucr  and  Solomon  (1973)  also 
measured,  by  means  of  a different  technique,  the  rate  for  HF(v  = 1)  and 
HF(v  = 2)  relaxation  by  F atoms  at  2545  K.  They  monitored  the  formation  and 
decay  of  HF(v  = 1)  and  HF(v  = 2)  in  the  reaction  of  F with  HC1  behind  an  incident 
shock  wave.  Their  estimate  of  8 x 10 * ^ cc/mol-scc  for  the  deactivation  of 
HF’v  = 1)  is  slower  by  a factor  of  6 than  both  their  previously  measured  value 
[Solomon  et  al.  (1971)]  and  their  estimate  for  HF(V  s 2)  deactivation  by  F.  This 
slow  rate  coefficient  agrees  with  the  value  obtained  by  Blair,  Brcshears,  and 
Schott  (1973), with  laser -induced  fluorescence  behind  a reflected  shockwave 
in  mixtures  of  HF  and  F^.  In  these  experiments,  F^  was  completely  dissoci- 
ated before  the  vibrational  relaxation  process  was  initiated  by  laser  pulse,  so 
that  the  F-atom  concentration  was  well  characte rizeG,  in  contrast  to  the  previous 
experiments  in  which  the  rate  of  or  SF^  dissociation  was  required.  Note 
that  Blair,  Brcshears,  and  Schott  (1973)  require  the  HF-HF  relaxation  rate  in 
order  to  extract  an  HF-F  rate.  Their  HF-HF  rate  is,  as  was  noted  previously, 
faster  by  a factor  of  1 . 5 to  2 than  other  measurements  over  the  same 
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temperature  range,  probably  making  their  calculated  HF-F  rate  a little  too 
slow,  but  the  effect  would  not  account  for  all  of  the  discrepancy  between  the 
different  workers.  Of  the  four  measurements,  the  faster  two  were  obtained 
under  conditions  for  which  the  underpopulated,  first-vibrational  level  was 
relaxing  toward  an  equilibrium  concentration,  whereas  the  slower  rates  were 
measured  for  an  overpopulated,  HF(1)  concentration  relaxing  toward  a lower 
equilibrium  concentration.  Quigley  and  Wolga  (1974)  examined  the  deactivation 
of  HF(1)  by  F,  O,  and  H atoms  at  room  temperature  using  laser-induced  fluo- 
rescence to  produce  the  HF(1)  and  esr  to  measure  the  concentration  of  atomic 
species,  which  were  generated  by  discharge  through  the  parent  diatomic  mole- 
cule. They  reported  a rate  coefficient  of  1.7  X 10*  1 cc/mo)-sec  for  HF-F 
deactivation;  O atoms  were  over  an  order  of  magnitude  more  efficient  at 
relaxing  HF(1).  Their  data  require  a small  correction  from  O atoms  in  the  F 
experiments  because  their  purification  of  the  was  limited  to  cold  trapping, 
which  would  not  remove  the  approximately  1%  O^.  -f  the  traces  of  O atoms 
react  with  F^  to  produce  FO,  the  effects  on  relaxation  are  difficult  to  predict. 
Several  theoretical  studies  of  the  effects  of  F atoms  on  vibrationaliy  excited 

HF  have  been  carried  out.  Thompson  (1972a)  performed  trajectory  calculations 

1 2 

and  obtained  rates  that  extrapolate  to  8 X 10  cc/mol-scc  at  2500  K,  agreeing 

with  the  lower  of  the  two  sets  of  experimental  numbers;  however,  the  laser- 

induced  fluorescence  data  between  1500  and  2500  K have  a much  steeper  tem- 

12 

peraturc  dependence.  Wilkins  (1973c)  calculated  a lower  value  of  2 x 10  at 
2500  K with  a different  potential  energy  surface.  Thommarson  and  Bcrend  (1974) 
tried  to  take  into  account  the  effects  of  multiple -potential  energy  surfaces. 
However,  they  concluded  that  their  nonreactivc  surface,  FHF  ^TT,  was  too  repul- 
sive, leading  to  an  overestimation  of  the  contribution  to  total  vibrational 
relaxation  rate  coefficient.  Shin  (1972)  assumed  an  inverse -power  potential  and 
calculated  analytically  the  transition  probability  for  the  decxcitation  of  HF(v  = 1) 
by  F.  His  calculations  extrapolate  to  9 X 10*  * cc/mol-ccc  at  2500  K;  however, 
his  rate  coefficient  increases  below  800  K to  about  2 X 10*^  cc/mol-scc  at 
300  K.  Neither  Wilkins  (1973a),  Thompson  (1972a),  nor  Thommarson  and 
Berend  (1974)  calculated  such  an  increase.  For  laser  modeling  purposes,  we 
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1 3.  2 

have  ^noscn  = 10  cxp( -2,  700/RT)  for  deactivation  of  HF  by  F atoms, 

which  has  the  temperature  dependence  obtained  by  Wilkins  (197  3a)  between  300 

1 2 

and  1000  K but  has  a value  of  9 X 10  cc/mol-sec  at  2500  K,  a compromise 
between  theoretical  and  experimental  results  (Figs.  9 and  10). 

• Wilkins'  (1973c)  theoretical  calculations  show  that  multiquantum  transitions 
contribute  significantly  to  the  deactivation  of  the  vibrational  levels  higher  than  2. 
Both  Wilkins  (1973c)  and  Thompson  (1972a)  found  the  rate  coefficients  for 
single -quantum  deactivation  HF(v)  to  be  between  v and  v^  times  as  fast  as  the 
deactivation  of  HF(v  - 1).  Rate  coefficients  listed  in  Appendix  1 .re  pro- 
portional to  v.  In  most  cases,  HF  deactivation  by  F atoms  will  have  a minimal 
effect  on  laser  performance  compared  to  other  mechanisms.  The  rate  coeffi- 
cient for  HF(v  = 1)  deactivation  has  a steep  temperature  dependence  and 
decreases  to  approximately  10**’^  cc/mol-sec  at  room  temperature.  Both  HF 
and  H atoms  are  usually  more  important  dcactivators . If  F atoms  arc  found  to 
contribute  significantly  to  the  deactivation  in  specific  lasers*  it  may  be  desirable 
to  include  reactions  involving  multiquantum  transitions. 

4.  M = F2 

Solomon  et  al.  (1971)  tried  to  measure  the  relative  efficiencies  of  Ar  and 

F^  in  their  shock-tube  experiments,  but  they  were  unable  to  measure  cither 

since  both  molecules  were  found  to  be  much  less  efficient  than  HF.  However, 

Bott  and  Cohen  (1971b)  were  able  to  obtain  results  for  M = Ar.  In  this  study, 

F2  Ar 

their  value  is  adopted,and  arbitrarily  it  is  recommended  that  k^  ^ z ^ 

Fried,  Wilson,  and  Taylor  (1973)  obtained  an  upper  limit  at  350  K of  2 x 

9 

10  cc/mol-sec,  which  is  consistent  with  the  foregoing  assumption. 

5.  M = H 

The  deactivation  of  HF(DF)  by  H atomirhas  been  studied  by  means  of 
classical  trajectory  calculations  by  Wilkins  ( 1973b)and  Thompson!  1972b)  as  well 
as  in  the  laboratory.  The  deactivation  can  occur  by:  (1)  the  loss  of  vibrational 
energy  to  translation  or  rotation  (nonrcactivc),  (2)  an  exchange  reaction  in 
which  the  incident  H atom  abstracts  the  F-atom  from  the  excited  HF,  thereby 
converting  the  vibrational  energy  to  kinetic  energy  in  the  liberated  H atom 
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(reaction),  and  (3)  for  levels  of  > 3,  the  reaction  of  the  incident  H atom  with 
HF(v),  forming  H->  and  F.  Wilkins  (1973)  calculated  the  deactivation  rates  for 
the  third  case,  H + HF(v)  l^v')  + F,  v > 3,  with  the  same  London-Eyring- 
Polanyi-Sato  (LEPS)  surface  used  for  the  pumping  reaction,  F + H->.  The 
deactivation  results  have  been  presented  in  Section  II.  A.  2.  Laboratory  ex- 
periments cannot  distinguish  between  (1)  and  (2). 

Three  independent  experimental  studies  have  been  reported,  all  at  room 
temperature.  Kwok  and  Wilkins  (1974)  worked  in  a large-diameter,  medium- 
pressure  flow  tube  at  room  temperature.  An  rf  discharge  through  SF^  gener- 
ated F atoms.  Downstream  of  the  discharge,  was  added  to  the  flow  and 
reacted  with  the  F atoms  to  produce  excited  HF.  Vibrational  relaxation  of  HF 
occurred  by  the  concurrent  processes  of  collisional  deactivation  by  HF,  H^, 
and  H.  Spontaneous  emission  also  contributed  to  the  de-excitation.  The 
observed  decay  of  the  various  vibrational  levels  of  was  simulated  by  com- 
puter with  the  various  individual  rate  coefficients  being  adjusted  to  give  a good 
fit  to  the  observed  data.  Experiments  were  performed  for  a variety  of  reagent 
concentrations  to  separate  the  coliisional  deactivation  of  HF  and  H from  that  of 

H?.  Kwok  and  Wilkins  (1974)  first  reported  that  HF(v  = 1)  is  deactivated  by  H 

^ 12  3 

atoms  at  a rate  of  less  than  1.9  X 10  cm  /mol-sec,  which  can  be  compared 

12  3 

with  10  cm  /:noi-sec  for  the  self-deactivation  by  HF,  which  is  made  in  equal 

3 3 

quantities  as  the  H atoms.  Much  faster  rates  of  1 . 1 and  1.8  X 10  cm  /mol-sec 
were  deduced  for  the  H atom  deactivation  of  HF(v  = 2)  and  HF(v  = 3),  respec- 
tively, with  the  HF  self-deactivation  rates  assumed  to  be  linear  in  v.  A more 

12 

recent  analysis  of  the  system  suggests  values  of  0.4,  6,  and  14  X 10 

cc/mol-sec  for  the  three  vibrational  levels  if  the  HF-HF  VT  rates  are  linear 

1 2 

in  v,  and  0.4,  4,  and  8 x 10  if  they  are  quadratic.  As  will  be  shown  later, 
we  believe  the  latter  assumption  is  probably  the  more  accurate. 


M.  A.  Kwok  and  R.  L,  Wilkins,  The  Aerospace  Corporation,  private 
communication. 
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A second  experimental  study  was  carried  out  by  Heidner  and  Bott  In  a 
very  slow  flow  tube  in  which  H atoms  were  generated  upstream  of  the  observation 
port  by  a microwave  discharge  through  II^.  The  atom  concentration  was  meas- 
ured by  isothermal  calorimetry.  HF  was  mixed  into  the  flow,  and  a small 
amount  was  excited  to  the  first  vibrational  level  by  a laser  pulse  from  a TEA  HF 
chemical  laser.  Deactivation  rates  were  calculated  from  the  decay  times  of  the 
infrared  fluorescence  with  and  without  H atoms.  The  results  of  this  study  indi- 
cated  a rate  for  H + HF(1)  — H + HF(O)  of  (1.4  t 0.4)  X 10U  cm3/mol-sec . A 
similar  experiment  was  performed  by  Quigley  and  Wolga  (1974)  in  which  the  H 
atoms  were  measured  by  esr.  They  reported  a rate  coefficient  of  < 1 X 10*^ 

Jjr 

cc/mol-sec.  However,  Heidner  and  Bott  argued  that  their  data  analysis  was 
faulty.  The  experimental  results  are  summarized  in  Table  11. 

Wilkins  (1973b)  calculated  the  deactivation  rates  of  several  vibrational 
levels  of  HF  by  H atoms  by  Monte  Carlo-selected,  three-dimensional  classical 
trajectories  on  a semi-empirical  LEPS  potential  energy  surface.  Rate  coeffi- 
cients over  a broad  temperature  range  could  be  deduced.  The  initial  trajectory 
calculations  were  performed  wi*h  a 1.5  keal/mol  barrier  height  in  the  potential 
energy  surface;  the  results  are  given  in  Table  12  (selected  results  are  shown  in 
Fig.  11).  Because  they  yielded  H-atom  deactivation  rates  that  were  considerably 
faster  than  the  experiments  suggested,  the  effects  of  the  barrier  height  were 
examined  in  a series  of  trajectory  calculations  performed  with  heights  of  2.5  and 
3,5  kcal/mol.  The  trajectory  calculations  for  the  three  sets  of  energy  surfaces 
all  predict  that  chemical  exchange  provides  an  efficient  mechanism  for  relaxing 
vibrationally  excited  HF  by  H atoms  and  that  multiquantum  transitions  are 
important  deactivation  processes.  The  calculation  with  the  1. 5 kcal/mol 
barrier  gave  a value  of  2.  5 X 10  cm  /mol-scc  for  the  room-temperatur e 
relaxation  rate  of  HF(1)  by  H atoms,  whereas  that  with  the  3.5kcal/mol  for  the 
surface  of  Polanyi  and  Schreiber  as  quoted  by  Bender,  Garrison,  and  Schaefer 
(197  5)  gave  a value  of  3.  1 X 10*^.  An  intermediate  value  of  3.  1 X 10*  * 


* 

R.  F.  Heidner  and  J.  F.  Bott,  The  Aerospace  C* *rpor ation , to  be  published 
in  J . Chem.  Phys  . 


64 


I n. 


3 4 3 4 7 9 9 to 
I0S/  T 

ID 


o i ; ) « * i r • 9 io 

io*/ 1 

<b> 


yjo  no  14?  .?*  too 


o I ? 3 4 3 ( ? • f io 

♦O1/! 

(Cl 


0 » ; 3 4 * 4 7 • 9 10 

♦oVt 

14) 

,,  t 'I  •' 

,.*».**♦*.  w,  MW>** 
•*  . *».*'*  ♦ «*, 


\ 


i 3 « 


\ r • i • o 


« * * ' # > 


Fig.  11. 


Rate  Coefficients 

f°rk6(V,v')for 
HF  (v)  + II  — HF(v')  + 

H Versus  10^/T  as 
Computed  by  Wilkins. 


(a)  v = 6; 
(c)  v = 4; 
(c)  v = 2; 


(b)  v --  5; 
(d)  v - 3; 
(f)  v = 1 

t 
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cm  /mol-scc  was  obtained  with  a potential  energy  surface  with  a 2.5-ktal/mol 
barrier  height  and  agrees  reasonably  well  with  the  experiments . For  each 
barrier  height,  the  trajectory  calculations  predict  that  the  overall  rate  coeffi- 
cient for  the  deactivation  of  HF(v  = 2)  by  H atoms  is  approximately  five  times 
larger  than  the  value  obtained  for  the  deactivation  rate  of  HF(v  = 1)  by  H atoms, 
which  is  somewhat  consistent  with  the  flow-tube  experiments  in  which  a much 
faster  rate  is  required  for  the  HF(2)  deactivation  than  for  HF(1)  deactivation. 

Recently,  it  wan  suggested  by  Bender,  Garrison,  and  Schaefer 
(1975)  that  the  H + FH  potential  surface  energy  barrier  is  no  less  than  40 
keal/mol,  which,  if  true,  would  rule  out  the  possibility  of  F-atom  exchange 
as  a significant  contributor  to  the  vibrational  relaxation  of  Hr  by  H atoms, 
leaving  only  nonreactive  processes  as  the  responsible  mechanism.  This  con- 
clusion is  consistent  with  experimental  results  in  this  laboratory  on  the  rate 

if 

of  the  exchange  reaction 


D + HF (0)  H + DF(0) 


6.  M = H2 

The  vibrational  relaxation  of  HF  in  the  presence  of  was  studied 
by  Bott  and  Cohen  (1972,  1973c,  1974b)  behind  incident  shock  waves  in  the  tem- 
perature range  1350  to  4000  K and  by  the  combined  shock-tube  laser-induced** 
fluorescence  (STL1F)  technique  from  2^8  to  611  K,  Because  of  the  near  reson- 
ance in  vibrational  frequencies  between  HF  and  (Table  15),  VV  energy 
exchange  takes  place  readily;  the  effect  is  to  couple  the  two  VT  processes:  HF 
deactivation  by  H2  and  deactivation  by  HF. 

k * 

6(1,  0) 

HF(  1)  + H2(U)  * HF(0)  + H2(0) 


R.  F . Hcidner.  Ill,  and  J,  F.  Bott,  The  Aerospace  Corporation,  to  be 
publ  ished. 
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The  experimental  data  cannot*  therefore*  readily  be  treated  analytically.  The 
high-temperature  shock-tube  results  were  simulated  by  computer  experiments; 
reasonable  matches  between  experimental  and  computed  profiles  were  obtained 
when  the  efficiency  of  H,  as  the  collision  partner  was  taken  to  be  0.3  to  0.5 
times  that  of  HF  itself.  However*  this  result  was  considered  an  upper  limit; 
the  actual  value  could  have  been  smaller  since  the  calculated  profiles  are 
relatively  insensitive  to  the  H->  efficiency  if  it  is  very  small. 

From  the  lower  temperature  STL1F  data*  Bott  ( 1974b)  was  able  to  extract 

HF  H>  3 -i  -l 

a value  fork^j  ^ ^§2.66k^  ^ of  (1  Z 0.3)  X 10  sec  Torr  , or  approxi- 
mately 6.2  X 10  T cc/mol-sec*  wheic  the  two  k's  are*  respectively,  for 
deactivation  by  HF  and  HF  deactivation  bv  H^.  **  impossible  to  isolate  the 

contributions  of  the  two  separate  terms  to  the  overall  expression.  Therefore, 

Hi  _ [ 

all  that  c.  n be  said  is  that  k . , € is  no  greater  than  approximately  390  sec 

-1  7 -t/3 

Torr  or  2.3  X 10  T cc/mol-scc,  OnaPt  versus  T plot,  this  expression  is 

consistent  with  the  one  derived  from  the  high-temperature  data.  The  contri- 
H F 

bunion  of  k-/f  is  expected  to  be  small.  For  example*  Bott's  (1974b)  overall 

'l1* 10 

relaxation  rate  coeificient  at  500  K is  3.  t X 10  cc/mol-scc*  whereas  the 

B 

H^-ll-*  VT  rate  coeificient  (discussed  below)  at  this  temperature  is  about  4 X 10 

cc /mol -sec.  HF  would  not  be  expected  to  be  much  more  efficient  at  chaperoning 

vibrational  relaxation  than  H,  itself.  ^Therefore,  at  temperatures  below  about 

1000  K,  it  is  reasonable  to  Assume  that  k^*j  is  described  by  the  expression 

6 X 10  T cc /mol* sec.  At  higher  temperatures,  this  expression  may  no  longer 

be  valid.  For  example,  at  2000  K*  the  latter  expression  gives  1.2  X 10** 

cc/mol-sec  for  the  overall  relaxation  Mtr  coefficient,  whereas  the  H,-H7  self- 

11  /,  z 

relaxation  rate  coefficient  is  about  1.6  X 10  . Therefore*  in  this  temperature 

regime*  the  contributions  from  Reactions  (6)  and  (7)  may  be  comparable. 


Bott's  (1974)  room-tempcraturc  value  is  in  favorable  agreement  with  the 
limiting  values  obtained  by  Hancock  and  Green  (1972b)  and  by  Osgood;  however, 
Hinchcn  ( 1973a)  has  deduced  a much  faster  rate  coefficient  from  his  data. 

7 . M = DF,  Dz 

Because  of  the  frequent  presence  of  in  combustion-driven  HF-lascr 
systems  (and, therefore,  of  DF  in  the  lasing  region),  there  is  interest  in  the 
deactivation  rate  coefficients  of  HF  in  the  presence  of  these  two  deuterated 
species. 

In  shock-tube  lasci>induccd-fluorescence  experiments,  Bott  and  Cohen 
.(197  3c)  studied  HF  relaxation  in  the  presence  of  DF  in  the  temperature  range 
of  295  to  1000  K.  HF  relaxation  occurs  simultaneously  by  VV  and  VT  processes 
with  DF  as  chaperone.  Therefore,  their  results  yield  a sum  of  two  rate 
coefficients,  for  the  two  processes 


H*  1 1 ) + DF(Q) 


DF 

'6(1,0) 


HF(0)  4 DF(0) 


and 


8(1,  G;0,  1) 

HF(  1)  + DF(Q)  HF(0)  4 DF(1) 


(8) 


DF  1 ** 

The  value  of  j Q Q ^ was  about  1.4  x 10  cc/mol-sce  at  300  K.  de- 

creasing slightly  to  (7  il)  x 10^*  between  550  and  1000  K.  The  room-temperature 
value  is  about  15%  larger  than  that  previously  obtained  by  Hinchcn  (1973a). 


The  experimental  data  can  be  described  by  a two -term  expression. 

k6U  0>  + k8(  1 0 0 1)  = 101,,03T'2  + I0*‘*T.  To  separate  the  VV.  Reac- 
tior  (8).  and  VT  Reaction  (6)  contributions,  the  VT  deactivation  of  HF(1)  by 
DF  can  be  assumed  to  compare  in  rate  to  that  of  the  deactivation  of 


R.  M,  Osgood,  MIT,  private  communication. 


6S 


of  DF(1)  by  HF.  The  latter  has  been  measured;  in  this  case,  the  experimental 

rate  is  a direct  measure  of  the  VT  process  because  the  VV  process  is  too 

endothermic  to  be  significant.  The  VT  rate  would  then  be  kP,*T  m = 10^‘^T  * 
i 3 3 ut  l , u i 

+ 10“'  T * making  a VV  rate  coefficient  k^j  ^ ^ ^ that  varied  between  3 and 
8 X JO1 1 cc/mol-scc  in  the  temperature  range  of  300  to  1000  K,  or  approximately 
6 x I011  cc/mol-scc  throughout  that  range,  within  50%  (Fig.  12). 

Bott  and  Cohen  ( 1 9*7 3c ) also  studied  HF  relaxation  in  the  presence  of 
at  the  same  time  by  the  same  technique.  Again,  their  results  y'eld  a sum  of 
rate  coefficients  for  VV  and  VT  processes.  is  a considerably  less  efficient 
dcactivator  than  s DF,  the  total  deactivation  rate  coefficient  at  room  tempera- 
ture being  (6  i 2)x  1 0 ^ cc/mol-sec  and  increasing  with  temperature  to  approxi- 
mately 1.3  X 1011  at  1000  K.  The  temperature-dependent  rate  coefficient  is 

9 0 6 5 

described,  within  10%  by  the  expres  jion  k = 1.38  X 10  T * cc/mol-scc,  where 
both  VV  and  VT  contributions  arc  included. 

B.  VIBRATIONAL- VIBRATIONAL  ENERGY  TRANSFER 
INVOLVING  HF 

Because,  in  general,  the  rates  of  energy  transfer  increase  as  the  mag- 
nitude of  vibrational  energy  converted  to  translational  energy  decreases,  VV 
transfer  in  which  only  a smalt  surplus  of  energy  is  transformed  into  transla- 
tional energy  is  very  fast  (much  faster  than  VT  transfer).  HF  is  unusual 
since  the  VT,  room-temperature  relaxation  rate  coefficient  is  only  about 
an  order  of  magnitude  smaller  than  that  in  VV.  No  experimental  d_ta  for  HF - 
HF  transfer  were  available  until  very  recently.  Because  theoretical  calcula- 
tions indicate  that  VV  transitions  in  which  the  quantum  numbers  of  the  collid- 
ing molecules  change  by  Av  > 1 arc  one  or  more  orders  of  magnitude  slower 
than  those  in  which  Av  - 1 [Rapp  and  Kassal  (1969)1.  it  has  been  customary  to 
i jnsider  only  VV  energy  transfer  reactions  of  the  type  shown  in  Reaction  (9) 


HF(v)  + HF(v') 


°( v.  v+  t;v\ 


1 l 


HF(v  + 1)  + HF( v'  - 1) 


<Q) 


T.  K 


1000/T.  1 K 


Fig.  12.  Vibrational  Relaxation  Rate  Coefficients, 
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Experimental  Overall  Relaxation  Rate 
Coefficient  for  MF{i)  + DF(0)  Versus  I03/1 
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If  the  molecule  losing  energy  is  an  harmonic  oscillator  and  the  perturbing 
force  linear  with  distance,  the  rate  coefficients  are  related  by  the  simple 
relationship 


k9(v,  vH;v\  v’-l)  = (v  + 1)v'k{0,  I; I, 


0) 


(a) 


Because  HF  is  very  anharmonicr  the  actual  relationship  among  the  various  k's 
may  be  considerably  different,  and,  in  fact,  may  decrease  rather  than 
increase  with  v. 

Another  effect  of  anharmonicity  on  the  value  of  k has  been  the  subject  of 
much  experimental  and  theoretical  investigation.  In  an  anharmonic  oscillator, 
the  VV  energy  transfer  is  not  a resonant  process;  some  internal  energy  must 
be  converted  to  translational  energy  during  the  reaction.  The  energy  converted 
to  the  translational  energy  is  given  by 


AE  a (E  *4  - E ) + (E  , i - E ,) 

' v+1  v ' v -1  v 

For  an  anharmonic  oscillator  in  the  lowest  rotational  state,  the  energy  levels 
to  the  first  approximation,  i.e.,  for  a Morse  oscillator,  are  given  by 

Ev  = hc*e(v  + i)-  ha»cxe(v  + I f (b) 

Thus,  the  energy  to  be  transferred  to  translation  in  the  above  VV  reaction  is 

AE'  = 2hcu)  x (v'  - v - 1)  (cj 

. e e 

The  reaction  will  bo  resonant  only  for  the  trivial  condition  tr.at  vT  = v + 1,  in 
which  case  no  reaction  is  physically  distinguishable*  For  greater  accuracy 
(seldom  warranted  because  of  the  current  state  of  energy  transfer  theory), 
expression  (c)  for  AE  can  be  extended  to  include  the  second-order  anharmon- 
icity teem  contribution 
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AE"  = 3hcu>eye(v2  - v'2  + 2v  + 1)  (d) 

Because  cuv  = 0 . 0 1 1 x for  HF,  A E*1  / A E*  = - 0»  0 1 6 5 ( v + v*  + 1).  For  v - 9 

c e c e 

and  v1  = 8,  the  correction  amounts  to  about  30%,  In  most  cases  of  interest,  it  is 
considerably  less. 

In  the  preceding  discussion,  it  was  assumed  that  AJ  - 0 for  the  energy 
transfer  process,  When  J is  allowed  to  change  during  the  collision,  the  amount 
of  internal  energythat  must  be  converted  to  translational  energy  can  be  de- 
creased consid,  rably.  For  example,  near  300  K,  more  than  20%  of  the  HF 
molecules  are  in  the  J - 3 level.  For  the  reaction 

HF(v  =1,  J = 3)  + HF(v  HF(v  = 0,  J = 2)  + HF(v  = 2,  J = 4 ) 

£E  is  only  85  cm”*.  The  large  possible  effect  of  AJ  on  the  magnitude  of  aE 
and,  therefore,  on  the  rate  of  an  HF-HF  VV  energy-transfer  process  may  result 
in  a very  weak  dependence  of  the  rate  coefficients  on  the  values  of  v and  v‘, 
i,  c. , on  the  value  of  aE  as  computed  by  Eqs.  (6;  and  (7). 

As  is  true  for  VT  transfer,  further  corrections  to  account  more  accurately 
for  the  effects  of  anharmonicity  vajry  from  theory  to  theory,  not  only  in  magni- 
tude but  in  the  direction  of  the  effect, 

t.  HF-HF 

Several  groups  have  reported  experimental  results  for  the  deactivation 
of  HF  by  VV  energy  transfer  [Bott  ( 1972),  Q3good,  Sackctt,  and  Javan  ( 1972, 

1973,  1974),  Airey  and  Smith  (1972)].  In  all  the  studies  reported,  it  is  necessary 
to  know  rate  of  the  competing  VT  process  in  order  to  extract  a VV  rate 
coefficient.  In  general,  the  competition  is  between 

k-9<v,v\0.  1) 

HF(v)  + HF(0)  ======  HF(v')  + HF(1)  (VV) 
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and 


HF(v)  + HF(O) 


HF 

C6(v,  v') 


HF(v')  + HF(O) 


(VT) 


Usually,  it  is  assumed  that  v'  = v - 1.  For  even  the  simplest  case,  i.c.,  v = 2, 
the  competing  VT  process  has  not  been  measured  directly.  Therefore,  its  rate 
coefficient  must  be  assumed,  and,  in  the  studies  published  to  date,  it  has 


HF 


= vk 


HF 


always  been  assumed  that  k^v  v 0)* 

In  a laser-induced  fluorescence,  room-temperature  study,  Bott  (1972) 

obtained  an  overall  rate  coefficient  for  the  disappearance  of  HF(2).  By  assuming 

k!^  = 2k!^T  = 2.48  X 10*^  cc /mol- sec,  he  obtained  0.95  X 10*^  for 

0(2,  1 ) 0(1,0) 

k f,n  M,  i.c.,  2 — 1;  0 — t,  or.  for  the  exothermic  rate  coefficient, 
i , u,  i ) up 

k^#J  ~ « = 2.2  x 10  cc/mol-sec.  If  k,  scales  faster  than  linearly  with  v, 

9(  1. 2;1,  0)  6 

then  the  VV  rate  coefficient  will  be  proportionately  smaller.  In  particular,  if 

,,  = 6k^  A.  * 7.5  x 1012,  then  kQ..  = 1 . 05  X 1013.  Note  that  the 

0(2,  1)  ^^-,0  ( 1 , 0)  ^2  9(1,2, 1,0) 

largest  n can  be  1.2  x 10  cc/mol-sec,  according  to  Bott.  It  is  apparent 

o(£,  i ) UP 

that  k^  2 i qj  is  only  weakly  dependent  on  k^^ 

Osgood,  Sackett,  and  Javan  (1974)  obtained  measurements  for  the  VV 
exchange  of  HF(v)  with  HF(0)  for  2 5 v < 4 in  experiments  in  which  one  to  three 
vibrational  levels  were  simultaneously  pumped  in  fixed  ratios  by  the  absorption 
of  laser  radiation.  The  exchange  rates  were  deduced  from  the  rise  times  of  the 
fluorescence  from  the  vibrational  level  above  the  highest  level  being  pumped. 

The  results  are  given  in  Table  13.  For  the  harmonic  oscillator  type  of  depend- 
ence previously  defined,  the  last  column  in  Table  13  should  read  1,  1.5,  and  2. 
Note  that  the  exothermic  rate  coefficients  for  v = J and  v - 4 (Column  4) 
exceed  the  gas  kinetic  collision  frequency. 

In  an  earlier  experiment,  Aircy  and  Smith  (197  2)  measured  the  relaxation 

of  upper  vibrational  levels  v = 2,  3,  4,  and  5 in  a steady-state  fluorescence 

quenching  experiment  and  deduced  values  for  krt/  , « that  increased  as 

M t>  e 9(v- 1 , v;  1 , 0) 

1:2.4:13.  5:44.  Again,  results  for  the  upper  rate  coefficients  exceeded  the 
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Table  13.  HF-HF  VV  Exchange  Rates  From 
Osgood,  Sackett  and  Javan  (1974) 

k 

• HF(v)  + HF(0)  -~-(v'  vr'1.;P.'  IWtv  - 1)  + HF(  1 ) 


V 

k 

exp 

( VV  + VT 
Deactivation)  |_ 

k-9(v.  v-l;l, 
HF 

0) 

HF 

6(1.0) 

k9(v-l.v;l,0) 

(Exothermic 

Direction) 

k9(v-l,v;l,0) 

6(v,  v - 1 ) v i 
assumed 

k9(1.2;l,0) 

2 

15 

12,  6 

30 

1 

3 

30 

26 

135 

4.6 

4 

26 

22 

260 

8.75 

Note:  All  k's  in  units  of  10*^  cc/mol-scc. 


gas  kinetic  collision  frequency.  This  experiment  is  somewhat  more  complicated 
than  that  of  Osgood,  Sackett,  and  Javan  (1974),  requiring  computer  modeling  as 
well  as  the  usual  assumption®  concerning  the  scaling  of  the  VT  rate  coefficients. 

The  large  VV  rates  deduced  from  the  data  of  Osgood,  Sackett,  and  Javan 
(197  2)  with  the  assumption  of  VT  rate  coefficients  linear  with  v (Table  13)  are 

inconsistent  with  observations  in  HF  lasing  systems,  both  pulsed  and  continuous. 

$ 

Modeling  studies  undertaken  in  this  laboratory  to  match  any  of  the  pulsed  HF 
laser  studies  of  Suchard  (1973),  Chen  et  al.  (1974),  and  Nichols,  Wrolstad,  and 
McClure  (1974;  show  that  when  the  fast  VV  rate  coefficients  are  introduced  into 
the  rate  package,  the  proportion  of  energy  that  comes  from  the  upper  vibrational 
levels  is  far  too  high.  That  is,  the  effect  of  fast  VV  rates  is  to  pump-up  upper 
levels  more  than  would  be  expected  on  the  basis  of  the  initial  pumping  distribu- 
tion. Similar  effects  have  been  observed  in  the  case  of  continuous  HF  lasers. ** 
The  conclusion  is  that  the  VT  rate  coefficients  must  scale-up  faster  than 
linearly  with  v,  and  the  VV  rate  coefficients  are  not  as  large  as  suggested  by 
Osgood  et  al.  (1974). 

*N.  Cohen  and  J.  F.  Bott,  The  Aerospace  Corporation,  to  be  published  in 
Appl . Opt. 

Emanuel,  TRW  Systems,  private  communication, 
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In  order  to  match  the  relative  energy  output  of  the  two  lowest  vibrational 
bands,  i.e.t  P-  (/P(  n,  from  all  three  of  the  pulsed-laser  experiments,  it  is 
necessary  to  have  ^ - ^^6(1  0)*  ^hc  behavior  of  the  VT  rate  coefficients 

for  higher  v levels  is  not  obvious  because  of  disagreement  among  the  experi- 
mental results.  At  present,  the  best  assumption  is  that  the  VV  rate  coefficients 

1 ^ hf 

are  constant  (**  10  ) with  v for  the  series  krt,  . - which  leads  to  k, 

9(v,  v-  1 ;0,  1)  6 

that  scale  as  1:8:28:25.  Experimental  results  ;.re  less  sensitive  to  other  VV 

coefficients,  i.  c.,  those  of  the  form  kn.  4 , for  v' >0.  Therefore, 

9(v,  v-  1 ;v\  v,J- 1 ) 

in  this  report  all  these  rate  coefficients  are  assumed  equal  to  those  for  which 
v'  = 0. 

Flow-tube  studies  in  a nonlasing  HF  medium  conducted  in  this  laboratory. 

❖ 

although  not  yet  completely  analyzed,  seem  to  point  to  the  same  conclusion. 

* 

Kwok  and  Wilkins  studied  the  vibrational  relaxation  of  upper  HF  levels  in  a 
medium-pressure  (1  Torr)  large-diameter  (10  cm)  flow  tube  in  which  HF(v) 
was  generated  by  the  reaction  of  F (obtained  by  passing  SF^  through  k dis- 
charge) with  H^*  Deactivation  processes  were  followed  by  monitoring  the 

vibrational  chemiluminescence.  They  obtained  overall  rate  coefficients  in 
1 2 

units  of  10  ce/mol-sec  for  the  disappearance  of  the  indicated  vibrational 
level:  k(v  = l...  5)  = 1 . 0:9.  9 ± 3. 0:15.  7 ± 6.  0:16.  3 ± 6.  0:5.  2 ± 3. 0.  These 
results  agree  within  a factor  of  2 with  those  of  Osgood,  Sackett,  and  Javan 

(1974)  for  v = 4.  Kwok  and  Wilkins'  additional  result  for  v = 5 suggests  that 

H F H F 

k/,r  is  considerably  smaller  than  k,  . . If  the  VV  rate  coefficients  are 

6(5,4;  7 6(4,3) 

constant,  the  VT  rate  coefficients  would  now  scale  as  1:8:28:25:5. 

2.  HF-H2 

Processes  of  the  type  shown  in  Reaction  (10)  are  expected  to  be  fast 
because  of  the  accidental  occurrence  of  near -resonance  b<  the  two 

molecules 


HF( v)  4 H2(vl) 


•I0(v.v+n;v'.  v'-n>  HF<V  + n)  + H 


(10) 


Kwok  and  R.  L.  Wilkins,  The  Aerospace  Corporation,  private  communication. 
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In  Table  14,  the  amounts  of  internal  energy  to  be  converted  to  translational 
energy  in  the  above  process  are  listed  with  all  species  assumed  to  be  in  the 

J = 0 rotational  level  or  all  species  assumed  to  be  J -2  level. 

i 

Recently,  Bott  and  Cohen  (1 973c),  using  the  STLJF  technique,  obtained 
VV  transfer  rates  for  HF  with  several  diatomic  partners,  including  (Fig.  13).  They 
found,  for  v = 0,  v*  = 1,  that  the  transfer  rate  is  nearly  temperature  independent 

Table  14.  Resonance  Defect  for  VV  Energy  Transfer 
From  H2  to  HFa 


HF 

H2 

Av, 

-1 

cm 

K b 
cq 

0 - 1 

1 - 0 

199 

(187) 

0.  38 

1 - 2 

1 - 0 

372 

(359) 

0. 17 

2-3 

1 - 0 

539 

(525) 

0.075 

3-4 

1 - 0 

702 

(688) 

0.034 

0 - 1 

2 - 1 

-36 

(-23) 

1.19 

1 - 2 

2 - 1 

136 

(120) 

0.  52 

0-2 

2-0 

336 

(305) 

0.20 

0-3 

3-0 

409 

(367) 

0.  14 

1 - 3 

2-0 

67  5 

(646) 

0.038 

aCalculated  with  all  molecules  assumed  in  the  J = 0 rotational  level. 
Values  in  parentheses  calculated  with  all  molecules  assumed  in  the 
J = 2 level  (the  most  heavily  populated  at  'I'  = 300  K).  Data  from 
Mann  (1961)  and  Hers&berg  (I960).  Positive  Av  means  an  exothermic 
process . 

^For  all  molecules  in  J - 0 rotational  levels. 
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Fig.  13.  Values  of  kvv  + k^  Obtained  for  the  Relaxation  of  HF 
by  Several  Molecules  Versus  T.  (The  data  for  H2 
were  multiplied  by  exp(288/T)  to  obtain  k_vv  for 
comparison  with  the  other  exothermic  rates.) 
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1 1 

over  the  range  300  to  1000  K,  with  a value  of  3.  3 x 10  cc/mol-ecc,  or  about 
a factor  of  20  slower  than  HF-HF  VV  transfer.  As  the  two  processes  transfer 
approximately  the  same  amount  of  vibrational  energy  into  translational  modes, 
it  is  apparent  that  the  rate  of  transfer  does  not  depend  solely  on  the  degree  of 
energy  mismatch. 

❖ 

Kwok  and  Wilkins,  in  a large-diameter  flow  tube,  studied  the  energy 
transfer  in  an  HF-H2  system  and  attempted,  by  computer  modeling,  to 
extract  numbers  for  the  rate  coefficients  for  processes  of  the  type  shown  in 
Reaction  (10).  An  important  qualitative  conclusion  of  their  work  was  the  con- 
vincing evidence  attesting  to  the  occurence  of  multiquantum  process,  i.e., 

An  > 1.  Quantitatively,  their  work  suggests  that  the  rate  coefficients  for  the 
procc*  Jes  written  in  the  exothermic  direction  increase  slightly  with  increasing 
v.  They  considered  only  processes  with  v*  = 0,  since,  in  their  studies,  the 
concentration  of  vibrationally  excited  was  too  low  to  contribute  measurably 
to  the  kinetics.  Their  results  arc  summarized  in  Table  15.  The  experimental 
observation  is  the  disappearance  of  a given  HF(v)  level,  which  is  a measure  of 
the  rate  coefficient  sum  given  in  the  penultimate  column.  The  breakdown  of 
this  sum  into  the  individual  contributions  indicr*td  by  varying  n is  a refined 
perturbation  on  the  total  disappearance  rate;  therefore,  ihc  values  for 
k4/w  . . are  subject  to  considerable  uncertainly,  c.g.,  on  the  order 

of  50%.  Because  of  this  uncertainty,  it  is  sufficiently  accurate  to  represent  the 
data  by  the  simplified  relationships  k10(v.„,  V;y'+n,  v’>  * vkI0(0,  ,;10)  = v X 
0.8  x 101  . The  temperature  dependence  is  assumed  to  be  negligible. 

3.  HF-N2 

Shock-tube  studies  of  the  relaxation  of  HF  in  the  presence  of  were  con- 
ducted both  by  Bott  and  Cohen  (1971a)  and  by  Blaucr,  Solomon,  and  Owens  (1972). 
The  former  observed  a slight  degree  of  non-Landau-Telle  r-type  behavior;  Pt 
increased  from  about  1 pscc/atm  at  3000  K to  about  3psec/atm  at  1500  K. 

Blaucr,  Solomon,  and  Owens  (1972)  found  tt  mporature-independent  rate  of 


* 

M.  A.  Kwok  and  R L.  Wilkins,  The  Aerospace  Corporation,  private 
communication. 
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Table  15.  Rate  Coefficients  for  HF(v)  + H2<v')  HF(v  - n) 
+ HF(v'  + n)a 


V 

V1 

n 

^10(v,  v-n;v',  \ 

V k 

/'+n)  n 1 0(v,  v-mjv1,  v'+n) 

10(v-n,  v;v'+n.  v* 
(Exothc  rmic 
Process) 

1 

0 

1 

0.  31 

0.31 

0.8 

2 

0 

1 

0.20 

0.44 

1.2 

2 

0 

2 

0.  24 

1.2 

3 

0 

1 

0.  18 

2.4 

3 

0 

2 

0.  12 

0.63 

3.2 

3 

0 

3 

0.33 

2.4 

aAll  k' 

s in 

12 

units  of  10 

cc/mol-sec. 

1*9  ± 0.5  psec/atm  throughout  the  same  temperature  range,  which  agrees  with 
the  results  of  Bott  and  Cohen  (197  la)  at  2400  K and  is  within  a factor  of  1 . 5 of 
their  results  at  other  temperatures.  Room- temperature  studies  in  two  labora- 
tories by  Hancock  and  Green  (1972a,  1972b)  and  Bott  and  Cohen  (1973c)  yield  a 
value  of  Pt  of  approximately  9 to  10  pscc-atm,  indicating  that  there  probably  is 
some  temperature  dependence  in  Prat  higher  temperatures. 

Initially,  these  results  had  been  interpreted  as  a simple  VT  energy  trans- 
ier  from  to  HF.  Later,  however,  Sentman  and  Solomon  (197  3)  e-examined 
the  data  and  concluded  that  while  deactivation  of  by  HF  is  probably  very  fast, 
deactivation  of  HF  by  by  VT  processes  should  be  negligible  at  temperatures 
from  1 500  to  3000  K.  They  concluded  that  the  above  experimental  studies  really 
measured  the  rate  of  VV  transfer  for  the  reaction: 

IIF(0)  + N2(1)^HF(1)  + N2(0)  (-11) 
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The  reported  relaxation  rates  are  ascribed,  therefore,  to  VV  transfer  at  this 
time;  the  VT  contribution  is  assumed  to  be  slight.  When  the  experimental 
data  are  converted  to  units  of  cc/mol-sec,  they  are  fitted  throughout  the  tem- 

3 9 

perature  range  by  the  expression  kj^j  q.q  jj  = HT  + 2.3  x 10  , In  the 
absence  of  any  experimental  evidence,  it  is  assumed  that  k^  ^ is  linear  with  the 
vibrational  level  of  HF. 

c.  h2  VIBRATIONAL  RELAXATION 

Occurrence  of  Hr -H2  VV  energy-transfer  reactions  can  permit  the  buildup 
of  energy  in  vibrationally  excited  molecules.  These  species  can  relax  by  VT 
transfer  from  to  irious  chaperone  molecules. 

H,(l)  + M — H2(0)  + M (7) 

Kiefer  and  Lutz  (1966)  measured  the  vibrational  relaxation  of  in  a shock  tube 

at  temperatures  from  1100  to  2700  K and  obtained  a Napier  time  of  3.9  X 
-4  1/3 

10  exp(100/T  ) psec-atm.  Their  H^-Ar  VT  rate  was  about  4.  1 times 

slower;  their  rate,  when  extrapolated  to  room  temperature,  agrees  wit’* 

a Napier  time  of  1.06  x 10^  psec-atm  obtained  by  DeMart  i and  Ducuing  (1966). 

When  converted  to  a rate  constant  expressed  in  the  form  ATn,  this  becomes 
H2  -343  r 

ky  = 10  T ; the  value  of  k^  is  1/4  of  this.  Recently,  Audibert,  Joffrin, 
and  Ducuing  (1974)  extended  the  temperature  range  down  to  40  K;  the  rate 
coefficient  below  approximately  100  K is  nearly  ‘?mperature-indepcndent, 
but  to  40  K.  the  above  expression  fits  the  data  within  a factor  of  approxi- 
mately 2. 

Heidner  and  Kasper  (1972)  studied  the  room- tempc rature  deactivation  ol 
by  A atoms  and  obtained  k^  = 10  * ^ cc/mol-sec,  making  H atoms  about 

1000  times  more  efficient  than  Itself.  Kapralova  et  al,  (1963),  in  experi- 
ments designed  to  elucidate  the  chain-branching  character  of  H--F-  mixtures 

M £ L 

(discussed  below),  obtained  deactivation  rates  for  k*  with  M = Ar,  He,  H~,  N?, 

H2  Ar  ' 6 c 

and  CF^.  Their  k^  /k  ratio  is  1.9,  somewhat  smaller  than  the  more 
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direct  determination  of  DeMartmi  and  Ducuing  (1966). 
efficiencies  were  all  within  a factor  of  2 of  that  of  H^. 
other  species  are  aseurruKi  to  have  the  same  efficiency 
of  dependence  on  the  vibrational  level  of  the  excited 
assumed. 


SI 


The  other  species 
In  this  review,  all 
as  Ar.  The  HO  type 
, i . e . , k oc  v#  is  also 


V.  CHAIN  BRANCHING 


The  possibility  of  the  occurrence  of  energy  chains  in  the  H^F^  explosion 
requires  brief  discussion.  Because  the  reaction  between  H atoms  and  F,  can 
produce  HF  in  highly  excited  vibrational  states,  it  is  possible,  in  principle,  that 
rw action  (12)  might  occur, 


Such  a reaction  would  be  exothermic  if  the  excited  HF  mo'cculo  can  transfer 
four  or  more  vibrational  quanta  during  the  encounter.  In  1963*  Vedeneev, 
Chaikin,  and  Shilov  (1963)  suggested  the  possiblity  of  the  occurrence  of  such 
reactions  to  explain  the  observed  upper  explosion  limit  in  m^x^urcs* 

Shortly  thereafter,  Kapralova  ct  al.  (1963)  claimed  to  have  obtained  proof  of 
chain  branching  in  the  H-*-F^  system  as  well  as  in  other  F,  systems. 

The  idea  of  energy  chains  is  quite  old  in  kinetics.  Such  processes  often 
have  been  used  to  explain  CO  oxidation,  H^  oxidation  sensitised  by  NO^  or  CINO, 
decomposition,  and  H%-Cl^  reaction.  As  pointed  out  by  Dainton  ( 1956 ),  in  no 
case  has  the  evidence  for  the  existence  of  such  chains  been  unequivocal.  (On 
the  other  hand,  ordinary  chain-branching  processes  arc  wcli-acccptcd  and 
occur  in  a wide  variety  of  systems.) 

Later,  Kapralova,  Margolina,  and  Chaikin  (1969)  favored  a more  elaborate 
branching  scheme,  i.c.. 


HF(v  * 4)  + F2  = HF(v*)  + 2F 


(12) 


HF*  + H2  - HF  + H-, 


(l  2a) 


H.>  + F,  a HF  ♦ H ♦ F 

A- 


(12b) 


S3 


because  their  experimental  results  did  not  support  the  one- step  energy  chain 
shown  in  Reaction  (12)*  For  the  branching  step,  they  obtained  an  approximate 
activation  energy  of  2 to  3 kcal/mol.  Alternative  explanations  for  the  effects 
of  inert  gases  that  did  not  require  postulating  energy  chains  were  offered  by 
Pavlov  (19&5)  and  by  Levy  and  Copeland  (1968)*  However,  Kapralova  and 
Chaikin  (1969)  presented  evidence  that  ruled  out  the  latter  explanations. 

Later,  Kapralova,  Margolina,  ;,nd  Chaikin  (197!)  attempted  to  answer  the 
question  of  which  of  the  two  processes.  Reaction  (12)  or  Reactions  (12a)  and 
(12b),  was  occurring.  They  argued  that  the  effect  of  pressure  on  the  upper 
explosion  limit  could  bo  interpreted  in  two  ways,  depending  upon  which  of  the 

two  mechanisms  was  responsible.  If  HF*  is  the  chain  carrier,  then  rieacti- 

❖ 

vation  is  due  mostly  to  radiative  decay;  however,  if  *5  the  carrier,  the 
fastest  deactivation  process  would  be  at  the  walls.  The  effect  of  varying 
vessel  diameter  suggested  the  latter  possibility*  As  a second  check  of  their 
argument,  they  added  a small  amount  of  HF  to  a mixture  and  found  that  it 
lowered  the  pressure  at  the  upper  limit  considerably.  T*  .s,  they  claimed, 

ijr 

would  happen  only  if  were  the  chain  carrier,  not  HF  . It  appears,  to  us, 

however,  that  the  effect  of  HF  could  be  the  same  regardless  which  species 

* t 

was  the  chain  carrier.  Furthermore,  since  it  is  likely  that  HF  and  H^,  if 
both  are  present,  will  be  in  rapid  equilibrium,  th*  n a process  that  removes 
either  of  these  will,  in  effect,  remove  both.  Therefore,  wall  effects  could 
hardly  be  used  with  confidence  to  distinguish  between  two  mechanisms. 

More  recent  work  of  Vedeneev,  Propoi,  and  Sarkisov  (1970)  on  the 
reaction  is  relevant  here*  By  the  use  of  arguments  si:  lilar  to  these 
of  Margolina  et  al.  (1971),  they  concluded  that,  in  this  system,  is  the 

important  chain  carrier,  and,  thus,  they  favored  the  analog  of  Reactions 
(12a)  and  ( 1 2b)  over  the  analog  of  Reaction  (12).  However,  they  continued 
their  discussion  by  pointing  out  cxj>erimental  observations  inconsistent  with 
Reactions  (12a)  and  ( l 2b)  and  concluded  that  the  branching  step  had  to  be 
Reaction  (l 2d') 

Dz(v  \)  ♦ F^  = F ♦ DF  + D (12'd) 
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where  the  D^v  = 3)  is  presumably  formed  by  a sequence  of  steps  such  as 


DF(v  > 1)  + D,  = DF  + 

D2(v  = 1) 

(12*') 

d2(v  = 1)  * D2^  - 1)  = D,(v 

= 2)  + D2(v  = 0) 

(12b') 

D2(v  = 2)  + D,(v  ^ !)=•- 

3)  + D2(v  = 0) 

(t2C) 

As  the  schemes  necessary  to  account  for  the  data  become  increasingly  involved, 
as  with  the  above  succession  of  mechanisms,  it  seems  that  the  probability  of 
their  being  correct  decreases.  Almost  as  an  afterthought,  Vedeneev,  Propoi, 
and  Sarkisov  (1970)  suggested  that  perhaps  there  is  a low-lying,  electronically 
excited  state  of  that  is  responsible  for  the  behavior  that  is  so  different  from 
the  other  halogens,.  The  implications  of  such  a process  arc  not  yet  clear,  but 
the  suggestion  is  worth  examination. 

An  alternative  to  the  mechanism  Reactions  (12a1)  through  ( 1 2d  * ) would  be  a 

$ 

sequence  in  which  several  quanta  were  transferred  concertcdly  from  HF  to  H^, 
but  there  is  no  evidence,  cither  experimental  or  theoretical,  that  such  a process 
is  likely.  Trajectory  calculations  that  could  provide  a better  assessment  of  the 

* 

probabilities  cf  Reactions  (12),  (12a),  and  (12a1)  arc  under  way  in  this  laboratory. 

Kapralova,  Margolina,  and  Chaikin  (1972)  subjected  some  of  their  explo- 
sion data  to  quantitative  analysis  in  order  to  extract  a numerical  value  for  the 
supposed  branching  coefficient  ^ j 2^*  When  some  of  the  rate  coefficients  values 

used  in  their  calculations  are  updated,  their  conclusions  lead  to  a value  of 
7 

k^b  = 1 . 2 X 10  exp(-4  100 /RT)  cc/mol-sec.  Because  of  the  simplifications  in 

their  kinetic  scheme,  necessitated  by  the  simplicity  of  the  analysis,  this  value 

4 • 

is  at  best  an  approximate  one.  This  expression  yields  a value  of  1.2  x 10 
cc/mol-sec  at  room  temperature,  whereas  the  deactivation  rate  of  H^(  1 ) by 
discussed  previously,  is  about  4.4  x 10^  at  room  temperature.  Thus,  in  a 
stoichiometric  mixture,  less  thar.  t/3000  of  the  excited  would  he  expected 

R.  L.  Wilkins,  The  Aerospace  Corporation,  private  communication . 
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to  undergo  branching.  However,  there  are  serious  difficulties  with  this  entire 
analysis.  For  example,  Kapralova,  Margolina,  and  Chaikin  (197  2)  deduced  a 
value  for  kj^g  j.j  qj  that  is  60  times  smaller  than  that  obtained  by  more  direct 
means . 

A more  extensive  computer  simulation  analysis  of  the  data  given  by 
Kapralova,  Margolina,  and  Chaikin  (1969)  was  carried  out  by  Sullivan,  Febor, 
and  Starner  (1975).  They  concluded  that  both  Reactions  (12)  and  (12b)  were 
necessary  to  explain  all  the  data.  They  deduced  limits  of  the  values  of  the 
corresponding  rate  coefficients  of  0.6  - 2 x 10^  for  anc*  ^ ~ 7 x * 0^  ^or 

kj^*  The  latter  is  about  a factor  o''  200  to  700  times  larger  than  that  deduced 
by  Kapralova,  Margolina,  and  Chaikin  (1972).  Even  if  th  sc  larger  values  arc 
correct,  they  are  unlikely  to  have  any  influence  on  most  chemical  laser  systems. 
Computer  experiments  carried  out  in  this  laboratory  indicate  that  the  occur- 
rence of  these  two  reactions  with  rate  coefficients  larger  than  the  upper 
limits  determined  by  Sullivan,  Feber,  and  Starner  (1975)  has  a negligible 
effect  on  H^-F^  pulscddaser  experiments 


N,  Colien,  The  Aerospace  Corporation,  unpublished  studies. 
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VI.  HYDROGEN  FLUORIDE  POLYMERS 


Because  of  the  great  strength  of  *he  hydrogen  bond  formed  between 
two  HF  molecules,  there  is  a strong  tendency  for  HF  to  form  not  only  dimers 
but  also  polymers  of  higher  order  in  the  gas  phase.  The  significance  of 
these  facts  lor  chemical  lasers  is  twofold:  (1)  The  attractive  force  that 

results  from  the  HF-HF  hy  Irogen  bond  is  probably  responsible  (at  least  in 
part)  for  the  increasing  rate  of  VT  deactivation  of  HF  by  itself  as  the  tem- 
peratuie  decreases:  and  (2)  The  presence  of  the  various  polymers  and  their 
possible  intluence  as  chaperones  promoting  vibrational  deactivation  must  be 
taken  into  account. 

The  first  of  these  two  issues  has  already  been  dealt  with  indirectly; 
the  empirical  expressions  for  HF  VT  self-deactivatior  include  the  contri- 
butions from  attractive  forces,  including  the  hydrogen  bond.  The  second 
matter  has  yet  to  be  treated. 

In  order  to  determine  the  contributions  cf  HF  polymers  to  any  relaxa- 
tion processes,  it  is  necessary  to  know  their  equilibrium  concentrauons . 

If  the  calculations  published  by  the  National  Bureau  of  Standards  in  Report 
No.  10 *904, Interim  Report  on  the  Thermodynamics  of  Chemical  Species 
Important  in  Aerospace  Technology  (1972)  are  accepted  as  reliable,  it  can 
easily  be  shown  that  under  most  conditions  of  interest,  dimers  and  hexamers 
are  the  polymeric  species  of  greatest  concentrations.  The  extent  Ci 
dimerization  and  hexamerization  as  functions  of  HF  pressure  and  tempera- 
ture are  shown  in  Figs.  14  and  IS.  For  example,  in  Fig.  14.  it  is  shown 
that  at  300  K and  an  HF  pressure  of  60  Torr,  the  fraction  of  H2  is  1%. 

The  same  fraction  of  H^F^  is  not  reached,  according  to  Fig.  15,  until 
[HF]  = 300  Torr  at  300  K,  but,  at  200  K.  both  H2F2  and  arc  1 To  of  the 

HF  pressure  at  0.3  Torr. 
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Fig.  15.  Hcxamorization  Calculated  Values  for 
[H^lF/-  )/[HFl  at  Various  Values  of 
[HF]  and  IoVt 
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If  the  contribution  of  polymers  to  vibrational  relaxation  is  significant, 

HF* 

then  the  apparent  value  of  k^j  qj  should  increase  with  HF  total  pressure. 

This  effect  has  been  observed  at  low  temperatures  by  Hancock  and  by 
Lucht  and  Cool.  In  Hancock’s  study,  the  experiments  were  performed 
with  DF  rather  than  HF,  but  since  the  thermochemistry  (including  the 
hydrogen  bond  strength)  of  DF  is  nearly  the  same  as  that  of  HF,  conclusions 
deduced  in  one  case  should  apply  reasonably  well  to  the  other.  Hancock's 
relaxation  data  included  a series  at  19?  K with  DF  in  the  pressure  range  of 
0.08  to  0.5  Tor r,  at  2 10  K,  0 to  0.8  Torr,  and,  at  300  K,  up  to  18Torr 
of  DF.  In  each  case,  he  observed  an  apparent  increase  in  k as  the  pressure 
of  DF  increased.  Calculations  in  this  laboratory  indicate  that  the  data 
could  be  explained  reasonably  well  if  it  is  assumed  that  the  dimer  has  a 
probability  of  deactivation  of  i/6  at  all  temperatures  under  consideration 
and  the  hexamer,  1/2.  It  then  follows  that  the  principal  contribution  at 
300  K is  from  dimer  (and  monomer);  at  200  K,  the  dimer  is  unimportant  and 
the  hexamer  (and  monomer)  are  responsible  for  the  observed  behavior. 

It  is,  perhaps,  still  too  early  to  try  to  fold  these  effects  into  HF 
chemical  laser  modeling.  The  uncertainties  in  the  kinetic  data  and,  more 
important,  in  the  thermochemical  data,  are  large,  and  refined  measurements 
need  to  be  made.  Another  important  question,  especially  in  cvv  lasers  with 
gasdynamic  effects,  is  whether  the  gas  has  time  to  reach  equilibrium  con- 
ditions with  respect  to  polymer  formation  after  nozzle  expansion  and  freezing. 
The  rate  of  polymer  formation  is,  therefore,  another  important  factor  that 
should  be  studied. 

Our  recommended  rate  package  at  this  time  does  not  include  values  for 
HF  deactivation  by  HF  polymers;  however,  we  feel  that  in  the  near  future 
these  processes  will  have  to  be  taken  into  account  if  accurate  laser  modeling 
is  desired. 


J.  K.  Hancock.  .Maval  Research  Laboratory,  private  communicat  ion. 
R.  Lucht,  Cornell  University,  private  communication. 
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N.  Cohen,  unpublished  studies. 
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RECOMMENDED  RATE  COEFFICIENTS  FOR  H7-F?  SYSTEM  (Continued) 
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APPENDIX  II 


** 


VIBRATIONAL  ENERGY  LEVELS  FOR  HF  AND  H2 


Vibrational 

Level 

V 

E 

V 

- E , cal/mol 
o 

HF 

h2 

1 

1 1. 327 

11,897 

2 

22, 161 

23, 121 

3 

32, 517 

33,687 

4 

42,  406 

43.  607 

5 

51,839 

52. 900 

6 

60,824 

61,  570 

7 

69,369 

69,640 

8 

77,477 

77. 120 

9 

85,  149 

84,020 

10 

92,  878 

90.360 

I 

\ 
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